{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Plot slice variable\n",
    "\n",
    "### David O. Lignell, May 21, 2020.\n",
    "\n",
    "Driver for the [get_var_slice.py](get_var_slice.py) function, which calls and extends the [slread.py](slread.py) script. The slread.py script reads a given slice file. This will be for a particular mesh, for a particular variable. The get_var_slice.py script provides a simple interface that will read all mesh files for a given variable for a given slice and combine them into a single variable for easy plotting or analysis. This is convenient for computing average or fluctuating fields, or otherwise processing the data contained in a slice. The interface requires minimal information from the user and reads the FDS input file to gather required information.\n",
    "\n",
    "This Jupyter notebook calls get_var_slice and plots the last frame along with the simple average of all frames.\n",
    "\n",
    "A complete case setup is available on [Github](https://github.com/BYUignite/fds_read_slice)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbYAAAGPCAYAAAAwfPUVAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzde5wsZX3v+8+vuntm1uIqWcpyy/Z2NGrUHBMwUQKKHlETdxKzTVSIIMaEHC8kJmabqJjAzo3NNgjH6Dkuj1vcJnhJNBEU5aIQLyARjAEMoGwRMbiAxX2tNZeurt/+46nqrqquvs309PT0fN+vV7+mu6u6qmYY5rt+z/PU85i7IyIiMiuijb4AERGRcVKwiYjITFGwiYjITFGwiYjITFGwiYjITFGwiYjITFGwiYjITFGwzTAz8yEfpwyx/8dzxz0jfS8xs8dWnHe7mT2c7nN+7v3Hl46ZmNm9ZnaxmT13hGuvupbs0TSz75vZ/2Nmh/b52Vi6n5vZCaVt9RF+dq8xsyelz/+gx7l+M91+VO69r5aOs2hm/2pmv2NmUenzTxpwDX/W6/sU2YrqG30Bsq5OKr0+FXgO8Bul96/KPf8S8OGKY32/4r1l4ATgv5Xe/yVgGxD3uK5PAJ8FasBTgTcCV5jZc9z9W6u8ljcDDwIHAP8XcBrw08AxPa7hWOBx6bFeA3wst61F98/uDcCRwG+W3v9a+n2sxg+Bt6fPd6TXcR7wY8CfVOz/MeDiivdvWOX5RWaSgm2Gufvf5F+b2YuAnym/X/LdAdvzLgZOpDvYTgQup3eofCt/DjO7CriQEB6/vcpr+ZS7706ffyCt6l5lZj/j7v9csf9rgH8DzgXeb2aPdPd7ADxMx1P+2b0U+D+rrsfMnjTkNZY9WPo5fAD4DvA7Znamuyel/b85ws9DZMtSU6SsxQXAT5rZ07M3zOww4KXA345wnC+lX58wxmv7Svr1/yhvMLN54FeBjwOfAhx41RjPvSruvghcCxxKqNpEZBUUbFK2YGY7Kh4LFfteRWjK+/Xce79GaIL8hxHOmYXPvWu4lrLHp1/vr9j2MuARwMfc/T7gEkIFNw0eDySEZtWy7T1+Ho3JXqLIdFOwSdlrgXsqHqdU7OuEqucEM7P0vV8HLnT3vX3Okf2BPtzMng+cn77/d2u4lsPSYz7OzF4HvCnd98sV+74GuM7db01ffwz42TU0Ka5WLRdOTzGzvwKeBVzk7isV+59J9c/j+RO7YpFNQH1sUvZZwgCGspt67H8B8EfA0WZ2B6Ff7ZcHnONd6SNzL/A77v7pNVzLt0uvrwd+w933599MR0r+AvDO3NufAfYTQvnMAdc+Tk8lBFPep4Hf6rH/LrrDH+Cb47wokc1OwSZl/+7ulw+7s7vfYGY3EkLhdkLT3xcGfOxDhEovBu4Evt+jQhnlWl6ZnnsHYUTkE4B9PfabB64rVWhfZf2DrbxG1PcJIRYRmmPfSbj+xR6f/+4o/21EtioFm4zDBcDvA7uBv3P35oD9b12HP9BfyUZFmtmFhIrtAjM7qjS6MOtLu6LqIGb2s+5+zYjnXkq/buuxfXtpv8y+/M/BzL5GqL7+FKi8J05EBlMfm4zDxwij+J5BCLkNlTY/ngH8FKFCA8DMHkdoKn0/YZBL+fEQxYEww7qLcE/fU3psfyphQMgdA677esJtBm82s8es4jpEBAWbjIG7fx/4XeB0OsPsN9rHgR8Q+v8yvw4Y8G53//vyA/g84d63kVoy0gr1MuAXyzOxmNlOQmh+1d0fGuJw/x2YI1TAIrIKaoqUsiebWdXQ93vd/fO9PuTu752Wa0mvJzaz84C/MrNfcPeLCcF2o7vf1uNjFxLuZ3sx1TN89PN2wu0P/2Jmu4DbgMcSZio5gCGbFt3922b2BeBUM/vz9HaEzE/3+Hl8392/OuL1iswsBZuUvTB9lF1HqGg207V8EPhj4O1mthv4CeAv++x/MWFAy2sYMdjc/UYz+9n0fK8jNM0+QLjd4L+6+7+OcLh3Az9PGASTH8xyQvoo+wRh8IuIABZmDxIREZkN6mMTEZGZMrFgM7M3mdn1ZvZQ+rjazF6W227pEiR3pkt4XJmfg1BERGQYk6zYfgj8IWEpkaMIE9/+o5n9ZLr9bcBbCf0KzwbuBi4zs4MmeI0iIrLJbWgfm5ndRxhNtoswA8Vfu/ufp9u2EcLtD9z9Axt2kSIisqlsSB+bmdXM7NXAgYQh0k8AdgKXZvukS3h8GTh6I65RREQ2p4kO9zezZwJXAwvAXuBX0rkGs/C6q/SRu4CeMzCY2amEVaGpUTtyOweP/6JFZFN6mPv3uPsjN/o6ps1LXrDd99xXXsN2NN+8fvkSd3/pmC5p7CZ9H9sthGU5DiEs9PgRMztutQdz912EZkwOtsP8Z6Pjx3GNIjIDLk8+eftGX8M02nNfwjWXHLGmYzQe/b92jOly1sVEmyLdfcXdb3X369z97cC3gN8jTJ4LcHjpI4fntomIiAy00fexRYQlRG4jBFi75EpXST6W0AcnIiIylIk1RZrZWcDnCDOcHwScCBwHvMzd3czOBd5hZjcD3yFMqLuXKZgtXkRENo9J9rHtJCzJsRN4kLBe1s+7+yXp9rMJ61m9D3gEcA3wYnd/eILXKCIy0xIS9idV6/rOjokFm7ufMmC7E9bQOmMClyMiIjNqo/vYRERExkrBJiIiM0XBJiIiM0ULjYqIbCEJsN/jjb6MdaWKTURExs7MHm1mHzGze8xsycz+zcyen9t+vpl56fH10jHmzey9ZrbHzPaZ2YVmNnDaFAWbiIiMlZkdCnwNMOBlwNMIS5LdXdr1cuDRuccvlLafC7wCOIEwYcfBwGfNrNbv/GqKFBGRcXsb8CN3Pzn33m0V+y27e+W0iWZ2CPB64HXufln63knA7cCLgEuqPgeq2EREZHQ7zOza3OPU0vaXA9eY2SfM7G4z+5aZvdnMrLTfMen275jZB83sUbltRwINisuZ3QHcxIDlzFSxiYhsIYk7+9e+wPQedz+qz/YnAm8E3gOcRVjV5b3ptr9Ov34B+DShkns88GfAl8zsSHdfJsxS1QL2lI59V7qtJwWbiIiMWwRcm67iAvAvZvZk4E2kwebuH8/tf4OZXUdoZnwZIfDWdHIREZFx+hHwb6X3bgIe2+sD7n4n8EPgyelbu4EaUF77beByZgo2EREZt68BTym99+OEiqySme0AHkMIRYDrgCbF5cyOIIyw7LucmYJNRETG7T3Ac8zsnWb2JDP7NeB3CKu3YGYHmtm7zey5ZvZ4MzsOuIhwO8A/ALj7g8CHgLPN7EVm9lPARwkrw1ze7+TqYxMR2UISjH3e9zawNXP3b5jZy4G/AN4F/CD9+v50lxbwTOBk4FBClXYF8MrSUmVvAWLgE4Rlzb4InOzurX7nV7CJiMjYufvnCItLV21bBF4yxDGWCTd2nzbKudUUKSIiM0XBJiIiM0XBJiIiM0V9bCIiW0iCsT9pbPRlrCtVbCIiMlMUbCIiMlMUbCIiMlMUbCIiMlM0eEREZAsJM4/MbfRlrCtVbCIiMlMUbCIiMlMUbCIiMlMUbCIiMlM0eEREZAtJPGJfMr/Rl7GuVLGJiMhMUbCJiMhMUbCJiMhMUbCJiMhM0eAREZEtJMFYci1bIyIismko2EREZKYo2EREZKYo2EREZKYo2EREtpAEY38yv6bHIGZ2hpl56bE7t/1AM3uvmf3QzBbN7BYz+73SMebTffaY2T4zu9DMjhjme9SoSBERWQ+3AMflXrdyz88BXgScBNwGPA/4oJntcfePpvucC/wycAJwb/qZz5rZke6eP1YXBZuIiKyH2N1399h2NPBRd78iff19M3s98LPAR83sEOD1wOvc/TIAMzsJuJ0QiJf0O7GaIkVEZD080czuNLPbzOzjZvbE3LavAr9oZv8RwMyOBp4FfCHdfiTQAC7NPuDudwA3EUKxLwWbiIiMaoeZXZt7nFrafg1wCvBS4LeAncBVZvZj6fbfAf4V+IGZNYF/Av7Q3T+bbt9JaLrcUzruXem2vibWFGlmbwf+M/AUYBn4OvB2d78xt8/5wGtLH73G3Z8zqesUEZllids4lq3Z4+5H9dro7p/PvzazrwPfI/x9Pwc4jVB5/RKhefF5wLvN7Pvu/gXWaJJ9bMcB7we+ARjwX4HLzewn3P2+3H6XEzoUMysTu0IRERk7d99rZt8Gnmxm24C/BH7N3S9Kd7nezJ4F/AGhOXI3UAN2APfkDnU48JVB55tYsLn7S/Kv047AB4GfAy7KbVru0+EoIiKbjJktAE8FriD0nTUojpIkfZ11j10HNIHjgQvSYxwBPA24atD5NnJU5EGEb+L+0vvHmNndwAOEdtd3uvvdk744ERFZHTN7N6Fg+QHwKOBdwAHAR9z9ITP7J+AsM9tLaIp8PnAy8DYAd3/QzD4EnJ3mQTbc/3pCq15fGxls5wHfAq7OvfcF4NOE+xoeD/wZ8KX0voXl8gHSDstTARbYvt7XKyKy6SVE7E/m1vs0RwAfo9OU+HXgOe5+e7r91YTmyL8FDiOE27uAv84d4y1ADHwC2AZ8ETh50D1ssEHBZmbnAMcAx+Qv0t0/ntvtBjO7jvANv4wQeAXuvgvYBXCwHebretEiIjIUd3/1gO27gdcN2GeZMMjktFHPP/Hh/mb2HsKd5C909+/129fd7wR+CDx5EtcmIiKb30QrNjM7D3gV8AJ3v3mI/XcAjwF+tN7XJiIis2FiFZuZvY9Qep4I3G9mO9PHgen2A83s3Wb2XDN7vJkdR+h8vBv4h0ldp4iIbG6TrNjemH79Yun9M4EzCEM9n0kYGXMooUq7Aniluz88oWsUEZlpiRv7W2u+QXuqTfI+NhuwfRF4Sb99REREBtFckSIiMlMUbCIiMlMUbCIiMlO00KiIyBaSYJOYeWRDqWITEZGZomATEZGZomATEZGZomATEZGZosEjIiJbSOLGYqux0ZexrlSxiYjITFGwiYjITFGwiYjITFGwiYjITNHgERGRLSQhYrGlmUdEREQ2DQWbiIjMFAWbiIjMFAWbiIjMFA0eERHZQtxhOZntP/2q2EREZKzM7Awz89Jjd277fzazS8zsnnTbcRXHmDez95rZHjPbZ2YXmtkRw5xfwSYiIuvhFuDRucczc9sOAK4Cfr/P588FXgGcABwLHAx81sxqg0482/WoiIhslNjdd1dtcPePApjZjqrtZnYI8Hrgde5+WfreScDtwIuAS/qdWBWbiIishyea2Z1mdpuZfdzMnjjCZ48EGsCl2RvufgdwE3D0oA+rYhMR2UISjKW1L1uzw8yuzb3e5e67cq+vAU4BbgYeBZwOXGVmT3f3e4c4/k6gBewpvX9Xuq0vBZuIiIxqj7sf1Wuju38+/9rMvg58D3gtcM46X5uaIkVEZH25+17g28CTh/zIbqAGlPvgDk+39aVgExGRdWVmC8BTgR8N+ZHrgCZwfO4YRwBPI4ym7EtNkSIiMlZm9m7gIuAHhD62dxGG+H8k3X4Y8Fjg0PQjTzKzB4Dd7r7b3R80sw8BZ5vZ3cC9hCbM64HLB51fwSYisoUkbiy11v1P/xHAxwhNifcAXwee4+63p9t/Cfhwbv8Ppl/PBM5In78FiIFPANuALwInu3tr0MkVbCIiMlbu/uoB288Hzh+wzzJwWvoYifrYRERkpijYRERkpijYRERkpqiPTURkC/HJDB7ZUKrYRERkpijYRERkpijYRERkpijYRERkpsx2D6KIiBQkGMvxbP/pV8UmIiIzRcEmIiIzRcEmIiIzRcEmIiIzZWLBZmZvN7NvmNlDZnaPmV1kZs8o7WNmdoaZ3Wlmi2Z2pZk9fVLXKCIy69yN5VZ9TY9pN8mK7Tjg/cDRwAsJ6+xcni44l3kb8FbCMgXPBu4GLjOzgyZ4nSIisolNLHrd/SX512Z2EvAg8HPARWZmhIXlznL3T6X7vJYQbicCH5jUtYqIyOa1kX1sB6Xnvz99/QRgJ3BptoO7LwJfJlR5IiIiA21ksJ0HfAu4On29M/16V2m/u3LbCszsVDO71syubbK8PlcpIiKbyob0AprZOcAxwDHu3lrtcdx9F7AL4GA7zMd0eSIiM8sdVlq1jb6MdTXxis3M3gOcALzQ3b+X27Q7/Xp46SOH57aJiIj0NdFgM7Pz6ITazaXNtxEC7Pjc/gvAscBVE7tIERHZ1CbWFGlm7wNOAl4O3G9mWb/ZXnff6+5uZucC7zCzm4HvAKcDe4ELJnWdIiKyuU2yj+2N6dcvlt4/EzgjfX42sA14H/AI4Brgxe7+8CQuUERENr+JNUW6u/V4nJHbx939DHd/tLsvuPvz3f3GSV2jiMisc4xmXFvTYxTprFNuZn/dY/sH0u1/UHp/3szea2Z7zGyfmV1oZkcMc07NFSkiIuvCzJ4DnApc32P7rwI/A9xZsflc4BWEcRnHAgcDnzWzgcmqYBMRkbEzs0OAvwV+g85EHPntjyPcz3wi0Kz47OuB/+Lul7n7NwljNH4SeNGgcyvYRERkPewC/t7dryhvMLM68DHgz9z9porPHgk0KM5EdQdwE0PMRDX90zSLiMi02WFm1+Ze70onzADAzH4LeBLwmh6fPxPY4+7/b4/tO4EWsKf0fs+ZqPIUbCIiW4g7Iw8AqbDH3Y+q2mBmTwH+gjCzVLNi+3HAKcCz1noRvagpUkRExum5wA7g22YWm1kMPB94Y/r8JcCjgR/ltj8O+G9m9sP0GLuBWnqcvKFmolLFJiIi4/SPwLWl9z4MfJdQyd1DGFSSdwmhz+2D6evrCANKjiedoCMd6v80hpiJSsEmIiJj4+4PAA/k3zOzfcB9ufuS7yptbwK73f2W9BgPmtmHgLPN7G7gXuAcwm0Dlw+6BgWbiIhMo7cAMfAJwoxUXwROHmZFGAWbiMgW4m60kskOr3D34wZsf3zFe8vAaeljJBo8IiIiM0XBJiIiM0XBJiIiM0XBJiIiM0WDR0REthAHWvFs1zSz/d2JiMiWo2ATEZGZomATEZGZomATEZGZosEjIiJbiRuJBo+IiIhsHgo2ERGZKQo2ERGZKQo2ERGZKRo8IiKylTh4a7Zrmtn+7kREZMtRsImIyExRsImIyExRH5uIyFbi4LFt9FWsK1VsIiIyUxRsIiIyU0ZqijSzeeA/ANuAe9z9nnW5KhERkVUaWLGZ2UFm9gYz+zLwIHArcCOw28x+YGYfNLNnr/eFioiIDKNvxWZmvw+8E/gecCHw58CdwCJwGPAM4FjgMjP7OnCau393Xa9YRERWz8G2+OCR5wDPd/dnu/ufuvsl7n6Du9/q7v/s7v/D3V8HHE4Ivuev+xWLiMhUM7M3mdn1ZvZQ+rjazF6W2+49Hu/L7WNmdoaZ3Wlmi2Z2pZk9fZjz963Y3P2VwxzE3ZeB9w+zr4iIzLwfAn8IfJdQQL0W+EczO9LdrwceXdr/KOAi4JO5994GvBU4BbgF+GNC6+BT3P3hfifXfWwiIjJW7v6Z0lvvNLM3AM8Frnf33fmNZvbLwHfc/Z/S1wa8BTjL3T+Vvvda4G7gROAD/c4/dLClIyLfCLwAeBSlZkx3/5lhjyUiIluDmdWAXwMOBK6q2H4g8GrgzNzbTwB2Apdmb7j7YjqI8WjGFWzAB4H/BHwG+DfAR/isiIhMBYO1Dx7ZYWbX5l7vcvddhbOYPRO4GlgA9gK/4u43VBzrRGAO+EjuvZ3p17tK+94FPGbQxY0SbL8E/HJWKoqIyJa1x92PGrDPLcCzgEOAXwU+YmbHufuNpf1+C/jMOO+LHmXmkbuBPeM6sYiIzC53X0lH0F/n7m8HvgX8Xn4fM3sWYeDIB0sfz/rgDi+9f3huW0+jBNs7gL8ws0eM8BkREREIeTNfeu9U4Dbg8tL7txEC7PjsDTNbINw33dVPVzZKU+SlwG8Dd5vZbqCZ3+juTxzhWCIiMqPM7Czgc8AdwEGEfrTjgPy9bNuBXwfOdvfCmA13dzM7F3iHmd0MfAc4ndBXd8Gg848SbP8T+AngXEIHngaPiIhsNg7WWveZR3YCf5N+fRC4Hvh5d78kt8+rgAOAD/c4xtmEeYnfBzwCuAZ48aB72GC0YDseeKG7XzPCZwrM7HnAHwBHEiZTfp27n5/bfj7hRr68a9z9Oas9p4iITJa7nzLEPh+md6iRVnFnpI+RjNLH9gNgedQTlBxImED5dwnzTVa5nHBXevb4hTWeU0REtpBRKrbfA842sze6+62rOZm7XwxcDO3qrMpy+a50ERGRYY0SbH9HGNFyi5ktA3F+o7sfPKZrOsbM7gYeAP4JeKe73z2mY4uIyIwbJdjevG5X0fEF4NOEoZ6PB/4M+FI6cWZXM6iZnUoYLsoC2ydweSIim9xkBo9sqKGDzd0/MnivtXH3j+de3mBm1wG3E4aIfrpi/13ALoCD7TCN0hQRkf6DR8zsoFEONur+g7j7nYTlD548zuOKiMjsGjQq8rtmdrqZHdFrBzOLzOznzewy4E3jvDgz20GY8PJH4zyuiIjMrkFNkccCfw58z8xuAK4F7gSWCDfM/QRhle1F4C/onu+rIF2e4Enpywh4bDpX2H3p4wzgU4Qgezzwl4Q5Kv9htG9LRES2qkEraH8XeKWZ/UfglYSg+xnC3eB7gH8h9HFd7O7JEOc7Crgi9/rM9PER4A3AM4GTgUMJ4XYF8Mph7jQXEZHhWDx4n81sqMEj7n4H8FfpY9Xc/Uqg33Ccl6zl+CIiIqPMPCIiIjL1FGwiIjJTFGwiIjJTRpl5RERENjlziOLZnnlEFZuIiMyUoYPNzL5kZu+qeP8RZval8V6WiIjI6ozSFHkc8Cwz+0ngZHfP1lObA54/7gsTERFZjVGbIl8EPAH4qpk9eh2uR0REZE1GDbYfAscAtwLXmtmR478kERFZNx5mHlnLY9qNEmwO4O5L7v4qwlRaVxKm2hIREZkKo/SxFcaHuvuZZvZt4PyxXpGIiMgajBJsLyDMwN/m7n9vZjcRJjcWERHZcKOsoP1PPd7/NvDtsV2RiIjIGugGbRGRrcQhaq3tMYiZvcnMrjezh9LH1Wb2stx2M7MzzOxOM1s0syvN7OmlYzzCzD5qZg+mj4+a2aHDfIsKNhERGbcfAn8I/DShq+pLwD+m90EDvA14K3Aa8GzCgtKXmdlBuWNckH7+penjp4GPDnNyzRUpIiJj5e6fKb31TjN7A/BcM7sBeAtwlrt/CsDMXksItxOBD5jZ0whhdoy7X53u89vAV8zsKe5+S7/zq2ITEZF1Y2Y1M3s1cCBwFWGSj53Apdk+6UxWXwaOTt96LrA33T/zNWBfbp+eVLGJiMiodpjZtbnXu9x9V34HM3smcDWwQAipX3H3G8wsC6a7Sse8C3hM+nwncI+7e7bR3d3M7k639aVgExHZQoyxzB6yx90H3eZ1C/As4BDgV4GPmNlxaz7zENQUKSIiY+fuK+5+q7tf5+5vB74F/B6wO93l8NJHDs9t2w080szaE4Okzx+V26cnBZuIiExCBMwDtxHC6fhsg5ktAMfS6VO7mtAn99zc558LHECx362SmiJFRGSszOws4HPAHcBBhNGOxwEvS/vKzgXeYWY3A98BTif0w10A4O43mdkXCCMkT00P+wHgs4NGRIKCTURExm8n8Dfp1weB64Gfd/dL0u1nA9uA9wGPAK4BXuzuD+eOcSLwXiD7zIXAm4c5uYJNRGQrcYjWeekZdz9lwHYHzkgfvfa5H3jNas6vPjYREZkpCjYREZkpCjYREZkpCjYREZkpGjwiIrKV+FhmHplqqthERGSmKNhERGSmKNhERGSmKNhERGSmaPCIiMhW4mCtjb6I9aWKTUREZoqCTUREZoqaImXiorm5VX0uWVnZ0GsY5/lFZP0o2GRiVhto5c+vJWDWcg3jOL+IrD8Fm6y7tQZar+ONEjDjvAYFnGxmBkQzPnhEwSbrpl+Y2FxjpGP5SrPn8fsFzLhDddTzi8jkKdhk7MYZaFWfK4dcVcAMCrTVXsdqA1ZEJkfBJmPVK1BWGyS9jtUvYAZ9dj3OnT+/Ak5kYynYZCxGqdJslc2DnguMfgGzXteQnT87Vr+AU7iJbJyJ3sdmZs8zswvN7N/NzM3slNJ2M7MzzOxOM1s0syvN7OmTvEYZ3bChZnNzIVDmGqt6tD9fcex+Rg61MZw/mptb1/49kVVziOK1PabdpCu2A4Ebgf+ZPsreBrwVOAW4Bfhj4DIze4q7Pzypi5ThjBJoQAgIgMYIv3bN3P9Fcw1YaWJzc0NXT32vY1Tp5yx96Ssrqt5EptBEg83dLwYuBjCz8/PbzMyAtwBnufun0vdeC9wNnAh8YJLXKr2NMjCjZ6g1SuHS7NGs2CME8+GSnTMfLpXXUHUdo2rGhYDNrkF9byLTY5r62J4A7AQuzd5w90Uz+zJwNAq2DTfK4IyhAy3T6/28ZjMcZ0C4DKzSyqE2zLl7yVWQw1RvoIATWW/TFGw70693ld6/C3hM1QfM7FTgVIAFtq/flclAI1VpjTpeD927FifDnaAZh89XhFvva+oTav3CrFc11+cays2jowxsEZHx2tSTILv7Lnc/yt2PajC/0Zcz09bcn9ZohOdpqHmjhjdqJNsa4XWfR9dx8q9LFWJXxZhXcS2VD6i+jvY+/a+h/DMp06AS2VAOFq/tMYiZvd3MvmFmD5nZPWZ2kZk9o7TPn5rZzWa2z8zuN7MvmtnRpX3mzey9ZrYn3e9CMzti0PmnKdh2p18PL71/eG6bTBGbaxSCpDLUGo12kHg9CkGWhtqgR6YQLNAdLP30Cliqw6v9qLieLITL31f+PPmfQ/7nU6Zwkxl3HPB+QjfSC4EYuNzMDsvtcwvwJuCZwDHAbcAXzCyfAecCrwBOAI4FDgY+a2Y1+pimpsjbCAF2PPANADNbIHwz/2UDr2vLq/oj3HdwRq8QyUKiXvyd9HoNi4ebvM6oF5sEs3NWNf0NCDWgEKCFa+rxfuE6qqTXoaZJ2crc/SX512Z2EvAg8HPARek+f1Pa5/eB1wPPAi4xs0PS169z98tyx7kdeBFwSa/zTzTYzOxA4Enpywh4rJk9C7jP3X9gZucC7zCzm1ffFtwAACAASURBVIHvAKcDe4ELJnmd0tvQow1zVVoWElmoJY3uhoLyOxa3CuFizc7rQrhBu6/LyAaR5K6r4nqya8lUhVg5fHtpB1zW7waFcIPetwXoVgDZQg4i/G9+f9VGM5sjjJd4CPhW+vaRQIPigMI7zOwmQiU4HcEGHAVckXt9Zvr4COHetbOBbcD7gEcA1wAv1j1s02GoUCsNEMk3K2ahljQikroRxd7zXFVt5JXhlp03fV5oDq2o1ICups7s2sqGCeD2tVRR9Saza4eZXZt7vcvdd/XZ/zxCYF2df9PM/hPwcWA78CPgeHfPBhDuBFrAntKx7qIz2LDSpO9ju5LOLUhV2x04I33IBus97+NooZaFRj7UvGFk4yGj2AshEjWT9uuoGfbKjmHNVvo6KlZL2bWsNEujIHOTJ+eqx857xdflMMte54O4Z7hVVJFV4Zanqk0mzZy+/6gc0h53P2qo85mdQ+hDO8bdy30OVxCaHncAvwV80sye6+4/WsvFTdPgEZki5VCrHASxllCrG61tEd4IX5N6eC+/b/k5VFdb5YqsfU255+3RlVCoIDPl8+SvN3/N2fvZZzvfZ2nUZGlQSfgZFkduZjSQRGaVmb2HMPDjhe7+vfJ2d9/n7re6+9fd/fVAE/jNdPNuoEYIvbyBAwqnafCITIHBs4r0GvlYPUikEFL5gKgVC/csUqqquKyCi6A9yKSySTK7riGaIPOBm1e+1ux5do1J7morK7c4qR6tWWp+LDdJ6uZtmTVmdh7wKuAF7n7zkB+LoH3v1nWEoDuedJxFOtT/acBV/Q6iYJO2QUvODOq/6hVq+UADSGpG0i5ahg+4/PasSTJcQ+7XuN3v1rsJshxq+TADugItH8LlcCs+S68tH26lQSW9miTbx1LAyQwws/cBJwEvB+43s6xPbK+77zWzgwlzA19E6Ft7JGHo/xHAJwHc/UEz+xBwtpndDdwLnANcD1ze7/wKNhlyqqzeoZYfadgv1LJAywIjyNr6Bwdc+xOlfrHC7CUDmiDLTaMwONDKIZyFW1KvDe53q3x/uPklFW6yXiYwQ/8b069fLL1/JmEMRQw8HfgN4McIofUN4Hnufn1u/7ek+36CMLDwi8DJFX11BQo2GV55UEZFM1+/UMtL2h9tT2dceh3kAy6p14gaEVEz6Q6S+hyWjY6sqCDD69rIVVoxiJ18uFnTC8fKX1O7osyN2Cw8H4LCTTYrd+85SDDdvh/4lSGOswyclj6GpmDb4kaZqR/o6j8a1PyYl4VEUi+HBQwTcNb09HOdIPF6DWu02kFi9L93btQqrf1+PbumYrgllPoDyz+/crNk9v6AJkkRWT0Fm4yu6ubrHjdel5sgs1Dz9Dcv6Rlo3QFXbJ4MAVqu3sqBBqtrdiwGGrnnqwy3TNdN3P3va1PVJjI6BZv0VBg0kvWtVcyKnx+Y0a9fLR9qWWB0wmq0gEsKr3LVW9xa1eCQXoFW7A/MrmcV4Va+5y6/TTdti4yVgm0LGzQKsqeKe9UGDRYph1pSg6gVnkcxQ1Zwvau3/PvjaHZsfyb3f4jF2ftrDLfcd6YZSWTSxnSD9lRTsMmaDBtqVZLSfdZZwGUBEsXeDr7sOOF/yO7qLT9Cca3NjuVQy1eX5SU7wjG6w201eoWbmiNFRqNgk4K+1Vqpb23QDPh55WqtsC1XvUHv5skoVzElDSPKMiALQULAwTBD94cPtDyvZ9cVAnbW/+Urshkp2GSwuUb1bBpr4DWw3J0oWbhBsXmyXL0FncCrDLiKG6tHCbT8tuzaMv3u/0nq1nN2ksIN5X2oSVJk7RRs0lflStRjUhVu0L96GybgejU1ts9TCLD+VVq5uix+rtjXFrV6V2+FKcAyQ8xEIiKjU7BJ28BBIz2spn8tUw436F29wXABFz43OMxguCqt/F5VX9uq9FggVVWbrCun08oxoxRsW9RI02j10OvetbJy/5qXQmOYcIPu5snwXlXAFT+Xv47OtZe3Fc9dvj7ovsayrFLMN0dmy+70o6pNZLy0bI2sWjbTfj9RHMLH4vR5qzsgqgIjKr03aG67qu39PtOv4iqfe1hRywuDSYYJNREZP1VssiadP95RGOre3mJETc8Nh3cSrHty46pQqwidfBC1Jx6ueC8vDDrJv+OFbflj1uLuKq6t1akgo5g0pL0d2lGTQv+aNV2jJUU2kIJNxsqaXmoGKA6yyMKtVzXVq5IqVEJx93uDPpMPuc79cBS2F8Kz61id6yuHWv582cTI46Z72USGp2CTVckPX8/P+5HUrR1uST0bKdgdblX6BdWg6qyqMzw/xL/4me6+uK6g63ON5VDraoKMPcxh2UywOEzQPOxwf5F15/QdwTsLFGwymmbcnkG/SnYfWa9wCzdZdw/sGNyHVvwfcZhRXeV9sqDrVckVVYRnrlKsCrVBTZAWJyMtWyMiq6Ngk+E1m5WTIEOnry1pRIWbpLPnWbgljYppsQY1KQ4TYhX/Ai2vAVcVdL3P3adi6xFqnfMMGDRSMRGyiIyPgk2Gk1vAE0pLsRCG/hf/oEeQhlqt6e0prqJWCJyoT1/UsM0kw/axdc/Qn52n87w7BLN5J/Pvda6vHGhZ82Nn32IzZPnnRTOuvIdNRNZOwSY9+cpKuJdtpdlZPTurNhr17j/WdPrbomaSrpeWTkocO5bexF2u6AYZZUBG/tiZWo/P5xdCLTdPtt/vuu3AC2EGnQotH2yVodaMK6s13cMmMl4KNhleVrX1CDdv1Ar3tlWGXBy1q7gqwwTdoKa+6oEkFX2CFSMjoTsIswDM96H1CrP286pQy6ifTTaQ+ezfjqJgk74qq7ZMGnRVlVtZPlaGnc5n1Buc88GSLTY66Fj5wOvug0u3DRlm0BktWhlq6lsTmQgFm7T5SnPwfJH5qq3R6Blu3qhVDnGv5Va4zhs0i8mow+Wr9q9aZqdWOm/+2nqFYbkyqzpvz1DLXqf9a1XNkJonUmRtFGwyugHh5vWo8Ae/V8j10mvftd4LZs1WzzXk2rPvDzFNWNW15F9XNj+KyMRorkgpqKoW2lVFfluvSmTQ8YdYnHTURUzzhmkWHfVcXq9VVpntc/YLtWYzfcSq1kQmRBWbjKYwQrJ35dbrBu6M12sDq6N84GTVVr+qLQuVYc5fPn6/6+w6T6+KclDzI2iIv2w8H22k8Wakik269K3aYKjKLV859QqCrBLqVw219x0QQu1ZPSrOv5rjDXtdQ/epQeHnpmpNZpmZPc/MLjSzfzczN7NTKvb5cTP7tJk9YGb7zeybZva03PZ5M3uvme0xs33p8Y4Y5vwKNqk08I/siOFWOHZFYAwbcuVAsjgpTlWVa/ZrbxvRMIEGFYG9hkpNoSYz5kDgRuB3gcXyRjN7AvA14DbghcAzgNOBvbndzgVeAZwAHAscDHzWzAb+D6qmSBlae+h/pl+zZNXnG53g6rtAacW2qhGK1myFgSpxkp477pw7nSUla5LMArH9dchFUgvXQBhcUm4S7QpW6Btqq7khWzP7y2bi7hcDFwOY2fkVu/w5cKm7vzX33veyJ2Z2CPB64HXufln63knA7cCLgEv6nV/BJj0NNfy/3+dzoZIPtXh7vec0V71EjYio2X9Fasv9Og8KtKQRreoa6vs74Qa5qi0L1lVQtSZbiZlFwC8CZ5nZF4Ajge8D73b3T6S7HQk0gEuzz7n7HWZ2E3A0CjZZizWHWy7U8oGSn86qSjl0PJtMeYiAy85b+Dri+SuviRDKWbjlj19o8hwh4BRqMnE+3Aw/A+wws2tzr3e5+64hP/soQlPlO4B3AX9EaI78WzPb6+6fA3YCLWBP6bN3pdv6UrDJSLqaI8saDWiEZW16hVprW9Q16XAv7aVmmkZU96EDrl+gJXUb+vx5dZJCuMHo99ZpXkiZEXvc/ahVfjbrA/iMu5+TPv+WmR0FvBn43FovTsEm6yofLK1tUTtU8rPmD2oSTOoQ1cN6aMMEXL4fryrQksbgc1b9izYmIqo7LCbE2+tEzYh8L3a7KVQ3Zov0sweIgX8rvX8T8Or0+W6gBuwA7sntczjwlUEnULDJQKM2R+artXy4ZMESb7OuhUYzXU2Q6X5JnAZbKeBaQBRH6QTF9faEy1XVWRZoYaHT/qt4V22vp6uARw0jnTWSqJm71y5rjhyyv03NkLIVufuKmX0DeEpp048TBocAXAc0geOBCwDSof5PA64adA4Fm4xHtlZbNhoxPwIyC5l2tZS+V/Hblw+8JF8OzUPUCqtdWxyWvomyr6Uqrqo6y58zqRePnV+axuPuQLM0o+KFEG5Jy8Lq4KQLq5IboamsEsHMDgSelL6MgMea2bOA+9z9B8DZwCfN7CvAl4AXEKq1lwO4+4Nm9iHgbDO7G7gXOAe4Hrh80PkVbDKSQv9auYpLh9pn/VtZtdYOmgbEC0ZrwXoGmOefZ9VaWqV5DK05sFa/kOuuzMph5jUKoZrQCa+kdIdM1KLzf8lSOGa8APWlsLxNWI4n6dwC0L7lIVe1zTU044hMDfPRV85YhaOAK3Kvz0wfHwFOcfd/NLNTCQNIzgO+C5ycDhzJvIXQZPkJYBvwxXSfgR3bCjZZnXyo5VbWbg+zz/rW0ibIeME6zZH1EFCdz3SeZ4HXDrU0aFrzIWQs7oRcUgtNkx4bLTohlw+z7BhZmCX19GupYrNcy2GUPrcYWul+1gIWAIz6Uvr5bREsJu2qLaxHl6TTizUrmyRtbk4DSGTmufuVQN+ObHc/Hzi/z/Zl4LT0MRIFmwzUt3+t3QTZ2afdt5ZvgkzDrbUQQi0fPFAMsk6l5u1tFqfNhPNpuOVCzuK0iktDLqvOqsIsf+zsuKGfztrhlu2TvY7i9L3lcI54Ifz/Wl8k9LfFnarNm63QHFkON1VtIhOjYJOhtZsh54qze+Rn+8gP7/eGEc9HxNs6oZaFS/YcusMGwGte+MocJC3DWtYOOasVQygLuarKLLwXju81LxzXWlY4ThZ0WYjmK8qY8D9NUk+bPFuG1Y2kEXWt7TZO0dycZh8RGZKCTdYmH2rpaMjCKMhGMWSS+RBqrfnwsXxVVgiz9H0WOn0BSTqww2KjFYeQA9phl1VdPUOs7nh6XMvOGxvE1h40YrVw3FbdQ8VW7wwosTQ4kzohzNLvqVO11dLmyNz0XhVVm5ojRdbXVAWbmZ0B/Enp7bvcfeCd5jIh5Wotfd4e4p+r1rImyKxfzWshFFrzxeABusLH6o7VEqJ6Qq0ewq0Vh/67JP3qrSgXOp2wy4IxH2JWC8fIjpU/ZhJHnWP1CbnWfLnLIO3Ty1VtFqfNkfQIt/yn5xoa8i+T5z6JwSMbaqqCLXULcFzu9fq178jq5Jsg8wNHGp3h/VXVWlapJXUnme+upLIAqqVhVosSGvUWjXr4FWjGtcLXVhICLh943oragQi0jxMuu1X4mh2rGddoJRGtOKK1XO8M+S+FHITRmF6HpNXpb4ua3VUb0B1uoL42kQmYxmCL3X33Rl+EFHVNo5UfLJJWa0BltRbPhybI0LeWhtp8UqioojSM5ubidpjN1VrM12Lm66HSWY47v67LrfB8pVUMu2ZcKwTYXK0TYvO1cJzseNkxl1t1Vloh4JbqSaeKq4cqzmLD605CRCvXLJk1ScbbQuDlq7b2z41cpbaGiZJFZHjTGGxPNLM7gWXgGuAd7v69AZ+REUT95nrsZy6t0Cr61bJmyHhbrXTPWmiCzKo1r4PPJ7CQYEBtPu6qzg6aW2a+HrNQyx6hwllqdcJ0qVWvfL4c19vBtVDrhEjxeTP32QZLrTpLrTrLcZ2H6/OFKq65VO/cwhAbrYUw+0h+xGSSVW5xp2rL60y11anaDM0bKbJepi3YrgFOAW4mzAB9OnCVmT3d3e8t75ze4HcqwALbJ3iZ0muiY6BdreWbIJNa2gQ5Fyq1qJbQWIgLgTZXa3Hw/FIh0A5t7M+ds/N0sVUM5+WKObrmo06Ybaut5J53gm2x1eCB5vYQcLUQig8tL7BSD+FWixJWVuokcUSyYHgr7S9Mv794PvTvZTeDW685KPN9kqraRNbVVAWbu38+/9rMvk5YfO61hOlUyvvvAnYBHGyHrXkdhq1g1dXaENrNkO1pszr3knkdkvkwqMPqPjDUDmkssq22wo819rG9tgzA/mwoJUBjX/vp/qT6e9oedVdE2bHan43CMRejOcLkBsD8Eg8tLxT2W6EeBpjUHK8bSQ2sDtFy/h68dB7J9IbtvNWs5F2mIf8yLraOt6ZMg6kKtjJ332tm3waevNHXMgvWFGr5Zshey9LUrfio5Yb51z0NBac+H/dsejx0bpH5KOawuX1sj1bY0XiYA6I0jOoPFy5pXzJfcaFF7c/2sC8Nv/1pFbdQq/PAyrbKcFvKVW2thXQgSa0T4kkDWO6sDr4e97Up3EQGm+pgM7MF4KkU5xyTCeu7/hqdOSGB9kz6nhven1VsWbVWq3eH2qFzS+2mx221JjsaD/PI+sNsj5ZZ6DGz8GG1vWv+3ranwZdVfYtpP17Wn7cc13mYEKCthZiVUtXWWgiznuTnprR6GCmZTQKdya/wrZGRIutnqoLNzN4NXAT8gNDH9i7gAMLEmbIGq63WqqbTKldrkJvBP98MmTVB5vrWsmrtgPkV5msxh6TNj/mmx6xKO6y2lwOiZQ6wUKFsj8YfBve0DmDBmtzXOhDoNE0+wHYOnYMHwgSRPEwYVGK1pF21ed3wVnHqrqRm1OhUbVmTZNeCpHMNWGXlpapNpL+pCjbgCOBjdBaX+zrwHHe/ve+nZHIKEx7XitVaqRmyPWdjqVrbvtAdauWmx8Nqe3lU7WG2R00OsBbbrXpQxnYb/Cu836sHa+x3h9q+dnBCp2kSsn63Trg14xor9aTd15bUPb0VIDcbScNDwANhvbbwrH1vWy6bwwwkqtxExm2qgs3dXz14LxnV2AaMNLqrt4w3+q9IPchiq1E52GMS9vn6Dqgp9LVly9qkVjv7iKo2WS1z725BmDHlwVsibf1m9c8qEOisNp3vX+uS3tjciiOacY3l9N4xKA7Vz/q6lrxRCJz9Ptyg1+3RcCGVHW9/0mifLzt/dg3LSe4ac/fKVS1Gmte9Cnitx54ish4UbDNuPYf3Q6cvaRTZLCHLcSc4FludQMlGO+5L5tmfNNiXrj46KNyqQq2qGTI7TnbccJ758DW9pWCxNVcYQJK/7rxs9e1sFYEs1EapYPsuC9THev+3FdmsFGxSqeuPbcWckEMfK1vnLFfpZFNhQWcE4mKr0Q6W/aWh/L3CrSq4+lVt+VALoTlXqNay68iqyKXS1F3ZBMzZygL95Cu3ws+sT5OuiKydgm2GTbJaazdHDui1zYIBqGyOzDcFZs2RWXNhv8ptmCbIcqhBsTrMzp/NalJuhswmXbZBTZG1zvbVVLSjUNUm0m2qBo/IbInizrpr1jKyOGoloZ9trtYKzZG1Ogu1JoutufZ0V/uSebZHyyFwIiAJw/33eY0DrMV+954jJatUhlrah5dVh/tb85X3sUF1M+Sowlpt1TOQaAkbmZgtMHhEwTajNupf8lEMSRyWdEkIfVCeDom3pYgEaNUTmlGNfczBPDyYzvCRn6g4GyHZnjkkDbf2c6CzolHc2ZbKmih7Bdq+ZJ4lb7A/meee+CD2t+a5t3kAi605Hmxu44GVhXBz9kpnUuTmUp1kuYal677Vlqy9+Gg0xN+Jfn9M1hJqGiEpUqRgk1WLmknfprbibFYGZMFWZ4VO5bZ9YYXlVp2D55fCLPtJnUMb+9mfzLGj8TD7kvnCDdv7WnMcYCvsp5HetN1ivzcLFdx+907TZY9A25fMs6d5EPuTOe5bOYDlJEyn9eDyQnspm/1Lc+212pLlGrYUEa1E1JaMqBWC3NKvUexEsaYsFdloCrYZtBHVWhR7mAQ4dqI4VGlJjXZFYzWoLaXhtuS0qJPECa10QdBmvcVKq8ZBc8vpMjINDmkssthq8GONfeyvzbUDbnu0zD5rhmouoR1w+xwOsFYh0LLmxl6Blp/dv1ylLa00CitsW2ztUGt/XwOqtUmtVKyqTaRDwSbrwvLNkWkLYyt935YjnHQarlZUqN4gDNZYjusszdc5dG4x7Xtbac8huT0KAbfkjULAQQg5oDLQsibHLNCyZselVr1dpe1bnqtcUdtiw5ajdpjVljvVWvb9ZqKWqjaRjaRgk+E04/YwdWu2+g75zxbfjGLaVVtEeuP2Ujp6kgjisIxNZfWW3tS81Kpz6NwSy0mdxVbMYqvBtlqzMP1WIeBSwwbaUqse1l/LVtAuVWksRe0+tWjZik2QucVG819FppqPZxmlfszsIOBPgV8hzPv7L8Dvuvs30u0G/AlhPc1HENbifJO7f3sc51ewzZhxNUP6SnOoG4ejZkKShlCoVCwNNSfB2lUbdAIuC7cwibAPVb1l80q2A66W3ntWmytUcAvWLDQ57k/m2NM8qB1oWT9atmL2UFVaKzQ9RnEItXYzZIv266xvLT9Pc76/rTBwpLk+ox/VHClT5P8HfpKwluYPgdcAl5vZT7j7vwNvA95KWFj6FuCPgcvM7Cnu/nD1IYen+9hkePk5Dput9mKFUexYM/3D3vJ2P1sUO7UlqC+HEKgthUfUCv1ttaVQCYUmPsNjo7VcD+GSROxfmmOl1Zl+K+t7K9zz1ireyJ3dbF3enoVafhh/OdSaS/WeoVZbMmrLne/BWuH7qi1R+H6zZsjs5wGzv6ijSJ6ZbQNeAfyRu1/p7re6+xnArcAb0mrtLcBZ7v4pd7+REIAHASeO4xoUbDNkIoNGmsX2tqiZtAdI5KuW/B97i9MQWAkjJWtLoY8qC4x2uC1FPcPtoeWFdrhBCKrsnrP9yVxh0dGsWsu25e9Ny45RFWreCue3pagdalkAZ6Fmcfg+6sudSq0damnud34OSffgkWZ3e+U471/TDdsyBepADVgqvb8IHAM8AdgJXJptcPdF4MvA0eO6AJFKvZojLU66+thC02P+X0oGhKbJKPb2si4sAMu5abZaYd+kDsk8sBThC0m7342FmP1Lc7Cw0r5ZOqu6ttVW2rOUbI9W2J/e1J3prtY6fWpZqK2s1Hv2p0WxtQM4ikMoR2nlmX3P+VDLqlVrdg/7DxXuZEZIikzADjO7Nvd6l7vvAnD3h83sauB0M7sR2A2cADyXULXtTD9zV+mYdwGPGcfFKdhmxKT/pV7oMyrdy9b5o56FW+75Uuh3i+dDQLQWwj/twvZOuFF3fAGaS3VYiGnGtfZK1ku1tGqL8qE2xwHRcmcWkXRqrmzex6xaW47r7YEi5VAjHc6f9ad1qspO02OhP60i1Lp+TnFr4rM8qK9N+nKvbDkY0R53P6rP9pOA/0HoX2sB3ySstXnkWk88DAWbrFmnuS3cgO0NK4Rb0sjCIIRbglGnM69ktsxNNgej142EJK3e0nBLLdfC0Hzml1hIA25/rcn22nL7PreqWfqzUHt4Zb5903WS3XQdWyHUspGP+abHbBRkVqVBMdQyWbVW2QxZ4isrmkZLZpK7/y/g+WZ2AHCwu//IzD4BfI9QwQEcDvwg97HDc9vWRH1sM2C9qzVfWYGVZvhXXjN8tTgJzWvpIJLsD3kUO7XFpPMHvuXUl7zQ71Zf8nbznuVGGYZh9CHcstnz8ysCNONa5aoAmQOKU50UZP1qmVZWqQGkzY89Rz7G3f1p9cVipVb4ntOfRVattZshcz8/1jnQ1Ncm08Dd96Wh9gjgJcBngNsIAXZ8tp+ZLQDHAleN47yq2Da59Q+13sP+y31GpQZJEqDW9NzaZEbUdJIG1JcgXrDCTdyteidQvEZYnLQeqqFWHFGbG9xHlU2cvL223O5/68dKQZq/Ny0/nL++5BQGiOT608L+nQEj4TilUBPZQszsJYQ/CTcDTwL+e/r8w+7uZnYu8A4zuxn4DnA6sBe4YBznV7DJaPJt8+kabfk+pHy4Rc2wbEsWcFHD0uVtQpNkNvUWgGc3PNdC1eR18JqH2wDqEdSHD4fQ11YdatkoyLws1MrVWlSaA7Lcl5YfJNIeGZprfizeu5ar1iZEfW2ygQ4B/hI4ArgP+BTwTnfPmirOBrYB76Nzg/aLx3EPGyjYNrWJNzetNCFfvTVjjHoIHjp/yL1ea0+QHP7Qh+3Zsyh2kvxCnQthXskoDv1t7cqtZni9enqqpVa9sBrAqJI4Hdqfvi7M/ZhrgqwthabUrGIL32fvQMvuWct+Fu0myLwJ9qsp3KSL+7pNEtA5hX8S+GSf7Q6ckT7GTsEmq9Nshim20nCDdL2xHtVb+53YSeoWAq6ZNkWm80pmVZvVKPxmhgmVI1r1MCPJcq3OfH1MlU82vD8uTpXVmVGkc1/eSIFWboJc5z8kItKhYNukJl2t+coKlp2zGYdmyCzc6D/3XFXAJUBUzwZjdFYDyKq2bA03z/WzDSsbFdkZ6l8c5t/KzdafDRrJP2or+RlFwgCY/MAQ6B9ohZ9H1gSZPc/9PCdFVZtsNQo2GU2+GS0Lt5QN+HUqDy7x2PtWbVlzJLGRUGOYmmd7NJ+7hy031L9Vb/evtZshq6q1XL9acTRncQj/UIEG3aG2QcP7FW6ylSjYNqGNGspdqNrKf6DTZkka9YGjALOAS+rWt2orrARAgtfDJMXNqMZKvRaWtqnVWajVWWw12B/Nt2fxv2/lgPZK2NlMI9nM/f2qtXITZH05obaYFIbwQ59Ag2KoVQTaJKs1ka1IwSYj6Qq3uUb1SL8BARcBUaO7aqsD7aMtFJsk831t+5eqw31brdm3CbK1XG/PMtKvWis0QZbuS4MBgQZDhdqkb85W1SbAuGYemWoKtk1mI6q1vkvY5MOt1DQJ9A24qBkVqrb8vW0DmyTTabaWayHAHljZBsBiK24vTfPg8kJ7NeyqJshB1VqhCXIxBMJQgZZ/TzOLiEycgk1GVqjaoLtyyw8syb3XdUN3BodZcgAADENJREFUvdau2kKoDR5IYlBoktxHeh3zS+2ZSNrrraX9alVNkNkkx8NUa4VKbdhAy34upZ/bRlPVJluBgk1WpWe4QXf11iPgrNEiaqY9btuingNJskmSLY5I5sPCpHHcIIkjWgudIFmOw20A5Tkhm0v1MCfkUtS1EnZtKczaH9ZW88pqrefAEChWqFMeaiJbhYJNhlLVHDkw3KA64NJws2aLKF19u9wkCZ0buLMVANqTJLdqJHMJCdAkTLfVjGtsX1hhuRX61PYvzbGyUm/3qVmfSY7zTZD1xe5qDfqMdsy/l/85lH5OVT9PEVkfCrZNZBontq0MN+iu3qDYPNmoh5uYGy1oRESxVzZJts9TJ6zltmS0FpwoW0kgNloLCStAK4moRUmxSsuthF1bsnaQ5Wfury3l5oJsrwDeuWetfdP5mEJto6k5cotzn/m+XwWbrFn2x3uk6g2gPtddtdGZkYQlb0+3BekCn2m4haVuIrzleGzEsZHMt2jVk85yNLmmxygOoVZb7swuki1HU5gLMu7uWyvMzJ9ZQ6CpWhNZXwo2GVrf0ZEMaJqEYvVGZyVuiztVWzYjSX0R4m0Wwi1dYbtFWmm1n5dW3o6NuB4mTs43PfZaCTubC7Ldr1Yxw0inWtvcVZrIVqJgk7EaOtyazXZzJJCr2mqFFarBwmCS9ONJbsJ8b4Xt0FmctD2UPx3On/Wn5VfC7gq1oZsgK+792aShpuZImWUKtk1iWvrXBlVtYZ8h+92y1QHSqi27abu9CkCzs7xNgnWm5Eqrr/ygkqzfzWve1Z82MNRyTZBA9aCRsormxGFCTc2QsuES3zT/AFstBZusi65wg4pbAtIJlEvL3UBEq2FErdAMGdUBvBBuCSGs8oNKvA5JnXZ/Wn6QSN9Q69EEWWiG7DOUP/t+RWQ6KNhkZMNUbT0Vwq1Z6nMLVVvSiELQpO9nTZJV4QadQSUJUEtHUuYHieRDrd3PVrFwaNd8kEOupTZsqE1btabmSJlVCrZNYFqaIUdVWbVBV+WWX88NaFdt5SbJQeFm9VDBtRcLrQy17tWw8+usQalaK193xfcoItNFwTbFpjXQRqnWeoZb1XGzgSTp6yQdUJI1SQ4KtwjaMyj3DrXO+fKhVq7WCprxzN73k/2OqXKTWaJgm1KzEGrTptwEOWnZz27amiRBAbeluE/l7+A4dS9uLBsqmpvbGqE2xP9Y+ebB1ShXa6tSXq1gDKb5HwfT+rsnMgpVbFNi2v+gTPMf47KsGXKa2Vxjav/VrOpNNjsF2xSY5lDb6ECzpkPdBu+4HtZ5McZpbpoEBZxsXmqK3EDT3OwI4wm1YQeObGUb/Y+HQab5d1Skiiq2DTLtfyzW/Y/tEMdP6oY3Rq/WPP2tzo+WTOjMYjK0qhXBS2xubixD/jdD9abKbTa4+8z/t5zKis3M3mhmt5nZkpldZ2bHbvQ1jdM0h5rNNTa8gkgaEUljdb+ayYB/qiV1I2lAUjOSYZo4G5P9WWz0z76faf69leliZm83s2+Y2UNmdo+ZXWRmzyjtY2Z2hpndaWaLZnalmT19HOefumAzs1cB5wF/AfwUcBXweTN77IZe2JhM8x+Hcf9RXU0zpKf3rg3eL4RYUgOvdV5D+OpDtkV4Y4iAa/Q/2LibW6c93Kb5d1imxnHA+4GjgRcS7jC93MwOy+3zNuCtwGnAs4G7gcvM7KC1nnzqgg34feB8d/+gu9/k7qcBPwLesMHXtSbT/gdhGv6YZjOPlCW1tMqqp8vU9AmifKB5LQTfoCoOWHWFuF6moXLuZ5p/l2XjuftL3P3D7n6ju98AnAQ8Evg5CNUa8BbgLHf/lLvfCLwWOAg4ca3nn6r/m81sDjgSuLS06VJC8m9K0/xHYL3+gI6rihmquTDVr0rzeicYe56rX7htUMgo3GRGHETIm/vT108AdpL7W+/ui8CXGcPf+mkbPLKDsBLJXaX37wJeVN7ZzE4FTk1fLl+efPLG9b28VVpa06d3AHvGcyEV1nZta7W+39vG0fc1HR630RcwjR7m/ksuTz65Y42HWTCza3Ovd7n7rj77nwd8C7g6fb0z/Vr1t/4xa7y2qQu2kaQ/yF0AZnatux+1wZc0drP6fcHsfm/6vmSauftLJ3k+MzsHOAY4xt0nMnXCVDVFEv412AIOL71/OLB78pcjIiKrZWbvAU4AXuju38ttyv6er8vf+qkKNndfAa4Dji9tOp4wOlJERDYBMzuPTqjdXNp8GyHAjs/tvwAcyxj+1k9jU+Q5wEfN7J+BrwH/N/AfgP9vwOf6te9uZrP6fcHsfm/6vmRLM7P3EUZCvhy438yyPrW97r7X3d3MzgXeYWY3A98BTgf2Ahf87/buLdSKKo7j+PeXUUmCWCaGXUQqwox6qB7EW0TSRej+4oOZZZElEUTQQ3QiLEPxqZeQ8ghJ0hWUvEAFBUVKFyojNCgzNZIKLMO04N/DmkO73fGcs8/MPnM5vw/Mw8yeWfNfzNnnz8ys/V+5zx8x8tN3DEbSMtJvHM4GdgEPR8T75UZlZmZDIelEieXJiOjJ9hHwBHAfMAHYATyQDf3Pd/4qJjYzM7PhqtQ7NjMzs7xqn9iaWFcyq58WbUvtRoVKmiNpk6QDWR8Wt33etVpx3TaEvvX2cw0/KincISu7xp9ZEWqd2BpeV3I36R1j33JpueEMyzjSO9KHgKP9fN61WnEjYLC+AbzNf6/hDSMTWi7zKLHGn1kRav2OTdIO4IuIWNqy7RvgtYh4rLzI8pHUA9weETMG27cuJB0BHoyI3mxdwEHguYhYkW0bS/on+UhEPF9WrJ1q71u2rReYGBELyoqrCJLGAYeBmyNic5OumzVXbe/YmlpXssW07FHPd5I2SppWdkAF62qtuIqYJemQpD2S1kqaVHZAwzCiNf7MilDbxMbAdSUn/3/3WtkBLAauA5aS+vOhpDPLDKpgA9WKq/v1A9gGLAKuIT22uwp4V9KppUbVuU5q/DXhulkDVPEH2qNeRGxtXc8GHXxLmtZhTSlBWUciYmPL6peSPgG+B24E3ignqs6UUePPrAh1vmMbNXUlI+II8BVwYdmxFKirteKqJiIOAvupyTUsq8afWRFqm9hGU13JrIbaxaQJV5uiq7XiqkbSRNJ0HJW/hmXW+DMrQt0fRQ63rmSlSVoNbAb2AZOAx4HTgfVlxtWpbETdBdnqScB5ki4Hfo2Ifd2sFddtA/UtW3qA10mJbCrwDGnk4JsjHWsnyq7xZ1aIiKj1AiwD9gLHSHdwc8qOqYA+bSQNqT4OHCD9g5xedlzD6Mc8IPpZerPPRUoAP5KmPH0PmFF23Hn7BowFtpMS2XHSu7Ve4Nyy4x5Cv/rrUwA9LfvU9rp5GR1LrX/HZmZm1q6279jMzMz648RmZmaN4sRmZmaN4sRmZmaN4sRmZmaN4sRmZmaN4sRmtSdpblZBf0zOdnZKuq2ouMysHE5s1gSrgBWRv1DvU8BKSf5emNWYv8BWa5JmkupovlJAc1tI849dX0BbZlYSJzarBEmLJP3SPl+ZpA2SNg1w6ELgnUiTXfYd0yNpl6Q7Je2V9IekdZJOkbRM0g/Zuda03p1ld3xbSAWAzaymnNisKl4l/T3e1LdB0njgFuCFAY6bDXzcz/apWVsLgFuBO4BNwJXAfOAeYHnWfqudwNzhdMDMqqHu1f2tISLiqKQNwBL+fay4EPgNeGuAQ88nFYxuNwa4KyIOA7skbSMlrCmRpjz6WtIHwNWkItN9DgJTJJ0cEX/n6pSZlcJ3bFYla4FrJZ2TrS8B1g+SYMaSKsy325cltT4/AXuypNa6bVLbcUdJ1etP6yhyM6sMJzarjIj4HPgUWCxpBnAF8OIgh/0MTOhn+1/tzZ9gW/t34Azgz0izlptZDflRpFXNWuBRYCLwQUTsHmT/z4DpBZ5/Bim5mllN+Y7NquZlYDJwPwMPGumzHZhV4PlnA9sKbM/MRpgTm1VKRPxOGjxyjKH9Nu0l4CJJl+Q9t6QpwExgXd62zKw8nkHbKkfSVmB/RCwd4v4rgbMi4u6c510FjI+Ie/O0Y2bl8js2qwxJE0iPAucDl3Vw6NPAckljcpbVOgSsznG8mVWA79isMiTtJY1KXBERz5YcjpnVlBObmZk1igePmJlZozixmZlZozixmZlZozixmZlZozixmZlZozixmZlZo/wDqT1GE9K+ZgsAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 576x432 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbYAAAGPCAYAAAAwfPUVAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzde7wsV13n/c+3uvc+JyEJECMXQQQFFBWezBAQMUBAkFFEGEVAJBAGzQiCqCDPA16ehAFGEQMZAeUwDgEUEBWBYAgQICICkSAIAcJFMFxyPSQETnL2pat/88eq6q6urr7t3fvW+/t+vfq1d1dVd63e55z9Pb+1Vq1SRGBmZrYosp1ugJmZ2Tw52MzMbKE42MzMbKE42MzMbKE42MzMbKE42MzMbKE42MzMbKE42HaIpDMkRfF4wIhjvlTsv3ibmzcVSadVPsOkx52nOP7XKu99cbHtSyPO/YDK686obD+j9p4dSV+X9FpJd5ih7U1tKR9HJX1K0m9KGvlvSNL3F8fn1XMX+x46w8/ujpJ+pfj+lBHn+ktJRyrP2w3vc6OkD0j66YbX/8qENjx01Oc0223aO90AYwV4AvBP1Y2S7gf8QLF/t/occHpt28uAK4E/rm2/Drhz8f0rgY82vN8ltecrwA9I+rGIqO/75WL/wRFtOwv492L/TwBPAh4k6Ucj4ubKcdO25SrgucX3J5H+zF4GfDfwuyPa8ETgMHBMcXz1Z3IZwz+7/wVcAfxJbfv1I95/Gu8FXg+I9PN/GvBOST8VEe9rOP4PgK80bL9sE20w21YOtp13AfCLkn4jItYr258AXA7kO9OsySLiGuAvq9skvRC4KiL+sn68pPLbD0XEm6c4xRWkz/8EKkEjaQl4DPAPwC+MeO27I6IMrP8t6Xrgt4FHAW+qHDdtW75d/UyS/pz05/NMSX8QEU1/Tr8MvAW4FSnkesEWEVcz/LP7Q+AbE352s/pCrd1vAz4FPAtoCrZ3RcSlGz2Z2W7grsid9ybgRODh5QZJLeBxwBubXqDkmZI+LWlF0rWS/kLSSbXjfk7S+UVX3KqkKyT9saSDtePOK97nDpLeJumIpOskvbRoy056I/C4Wjv+C+ln1vjzGeH9xde7zKNREbECfAw4HrhNfb+k+wJ3B95M+jO+l6Qfnce5NyMiPg18i9QbYLaQHGw77+ukbsgnVLY9lPTLctQv7j8DziFVMc8CDpEqmA/UQuspwCqpi+s3SL/cfws4r+E9M+BC4JvAc4B/BJ4NnLmBzzTJ8ZJOang0hegbgdsCP1nZ9gTgn4GvznDO8hf5NzfRlro7A0EKironAl8DPgS8m9Sd+MQZ2rslJN0aOAG4YcQht2z6eWxjE802zV2Ru8MbgXMk3SIibiJ1YV0SEf9e74KSdH/gvwNPjojXV7ZfSArIJ5GCDuCXa+NJr5b0ReCFkn4nIr5W2bcE/E1EvKB4/ueS/hV4KilI5+lQpY1V9yB17/VExFckfZT0M3mPpOOAnyOF7ji3LH4hl2NsfwAcBd65wba0Kr/gv4v0czkF+IeIOFp9oaQ2qeJ+Q6RVxtcl/R3wBEnPi+1defxg0e5yjO1FpP/E/O2I4y9q2ihpKSI6W9JCszlzsO0OfwP8KfBoSW8FHs3oCQmPBY4AF9b+J305cA3wYIpf1GWoFTP3jieF14dIv+T+M6miqHpN7fk/MTzBYR5eBFzcsH1UBfZG4EVKMxUfTRHCjO9WvLD2/LPAb0TENzbYlruSJsBUvZ0UcHU/Raq4q2N5bwJ+FXggqRreLk9lsI1rwIuBc0cc/+vAFxq279qxXrM6B9suEBHXS3o3qSrpAMcCfz3i8LsDx5FCrElvvKcY03kJcBppZl7VLWvP1yPiqtq2G4BbT2r/BlwWEY2VwQh/TZqB+EjSz+jdEfFNSeOC7TdIszZXSCH1tRGV0rRt+Rrw30j/KfgB4Pmk2ZFHG459Imlm6I2S7lps+wbp5/lEti7Ymj7f24FXAMvAfUjtXh5TNf6LJ4/YXudg2z3eSJqWfQJwUURcO+K4jDRO9PgR+28AkHRL4APATaTq70ukX8J3II2x1cdXu5to+5aKiGslvQ94JnA/pqsiP1aZFTkPN1cC8L2S/gn4JKn6+c3yoKKr9FGk/5x8seF9HiPpGRGxOuP5y8s+6v9BKR1L86UhX6+0+wJJh4FXSPrHiKh3y5otBAfb7vF20kSPnwCePOa4fwceBnw0Io6MOe7BpIriMRHRqxAkPWwObd0JbyQF8hHgHTvbFIiIz0h6HfA0SedERNl1+fOkkPl1oP6fkzuSKs9HAG+d8ZRXFF9/kNo1j4UfqhwzzqtJE4heLOkftnm8z2xbeFbkLlGMhz0NOBv4+zGH/jXpz+0P6jsktYpZb9AfE1Flf0a6lmsveivpZ/P02oSYnfTHpPG+6kSWJwJXRMSrIuJvqw/SuNZVbGx25L+QKvVfkXSguqNYFeQepOv6xiomgJwD3BP42Q20w2zXc8W2i0TEG6Y45oOSXgn8jqR7kaaSr5ImNzyGFHjnkabDfxN4naQ/BdaL/cdtTetncmoxc7DusxHxr00viIjvkFYT2fG2VNr0eUnvJIXN/yCF3EMYMYs0IqI4/kmSbhURTZcJjDrXqqTnAK8FLpX0ZtKqJvciTUr5EqkanMZrgRcAzwPOr+37aUk/1PCaSyKiqWvVbNdxsO1BEfGMYir+r5Fm9XVIEyTeQnEhcjEh5RGk5ZnOJnXh/R3pl+6ndqLdFb9ePOr+BBgbJruwLS8lTWp5Fml8s8VwWFS9nRREj6X5MoORIuI8SdeQrjP8HVKX55XF+5w9bVBGxFFJrwJ+X9KDql3VpMBr8jSaxwzNdh25i93MzBaJx9jMzGyhbFuwSfp1pVt9fLt4fKToKiv3S9JZkq5Uui3IxZJ+ZLvaZ2Zmi2E7K7avA/8vacWLU0hjQW8rJkBAuiXIs0nXKt2HNFX6vZKO38Y2mpnZHrejY2xKtxJ5Hmnw+0rgFRHxomLfMaRwe05EvHrHGmlmZnvKjoyxFddbPZ409fzDpDX/bge8pzymWFj2g8D9d6KNZma2N23rdH9J9wQ+Qlpx/QjwXyPi08WK9TC8/uE1pCWgRr3fmRS3VWnRuvexnDD/RpvZnvQdbjgcEd+90+3YbR7+4GPj8PWbW0HvXz+1+u6I+C9zatLcbfd1bJ8HTiYtwPsY0sXDp230zSKid8uRE3Ri/Fi2V1eLMrN5u6j7lmmWGNt3Dl/f5ZJ333FT77F0+3/f1ffo29auyIhYi4gvRcTHI+J5pEVkfwu4ujjktrWX3Layz8zMbKKdvo4tAw4AXyEFWK/kKu4E/QDSGJyZmdlUtq0rUtIfkhZp/RrpppdPIN0n7BHFGnovB54v6XLSjQ5/jzQO98btaqOZ2aILgtVY3+lmbKntHGO7HfCXxdcbSesV/nREvLvY/xLSvaZeSbq55SXATxWL35qZmU1l24ItIs6YsD9Iq7eftQ3NMTOzBbXTY2xmZmZz5dvWmJntI12Cm6Oz083YUq7YzMxsoTjYzMxsoTjYzMxsoTjYzMxsoXjyiJnZPtIlWInNLYK827liMzOzheJgMzOzheJgMzOzheIxNjOzfSRCrIR2uhlbyhWbmZktFAebmZktFAebmZktFI+xmZntI11gJRa7plnsT2dmZvuOg83MzBaKg83MzBaKg83MbB/pIm7uLm3qMYmk50n6mKRvS7pO0vmSfnTM8a+WFJKeU9t+QNKfSjos6SZJ75B0x0nnd7CZmdm8nQa8Crg/8BCgA1wk6cT6gZIeA9wXuLLhfV4O/ALwS8ADgBOAd0pqjTu5Z0WamdlcRcTDq88lnQ7cCPwEcH5l+/cB5wIPBd5Ve80tgacCT4mI91be54ri+HePOr8rNjMz22rHk/LmhnKDpDbwJuCFEfG5htfcG1gC3lNuiIivAZ8jVYIjuWIzM9tHuoiVmDxONsFJki6tPD8UEYfGHH8u8EngI5VtZwOHI+LPRrzmdkAOHK5tv6bYN5KDzczMZnU4Ik6Z5kBJ5wCnAqdGRF5sOw04Azh5KxrnrkgzM9sSkl5GmvjxkIj4cmXXacDtgaskdSR1gO8D/kjS14tjrgZawEm1t71tsW8kB5uZmc2dpHPph9rltd2vAu5FqtjKx5XAy4CfLI75OLAOPKzynncE7gF8eNy53RVpZraPpPuxbe2vfkmvBE4HHg3cIKkcEzsSEUci4lrg2tpr1oGrI+LzqZ1xo6S/AF4i6Vrgm8A5wKeAi8ad38FmZmbz9vTi6/tq288GzprhfX6TdA3cXwPHFO/3pHKsbhQHm5mZzVXE7Lfojog7N2xbBZ5ZPKbmMTYzM1sortjMzPaRdB3b8k43Y0u5YjMzs4XiYDMzs4XiYDMzs4XiMTYzs32ki7ipe2Cnm7GlXLGZmdlCcbCZmdlCcbCZmdlCcbCZmdlC8eQRM7N9pBtipbvpG43uaq7YzMxsoTjYzMxsoTjYzMxsoWxbsEl6nqSPSfq2pOsknS/pR2vHnCcpao+PblcbzcwWXSBWYmlTj91uOyu200i3A78/8BDSzeMuknRi7biLgNtXHj+zjW00M7M9bttmRUbEw6vPJZ0O3Aj8BHB+ZddqRFy9Xe0yM7PFspNjbMcX57+htv1USddK+oKk10i6zQ60zczM9qidvI7tXOCTwEcq2y4E3gp8Bbgz8ELg/ZLuXdwifICkM4EzAQ5y7Fa318xsz+siVhf8OrYdCTZJ5wCnAqdGRF5uj4g3Vw77tKSPA1cAjyAF3oCIOAQcAjhBJ8aWNtrMzPaEbe+KlPQy4JeAh0TEl8cdGxFXAl8H7rYdbTMzs71vWys2SecCjwMeHBGXT3H8ScAdgKu2um1mZrYYtvM6tlcCTwGeANwg6XbF47hi/3GSXirpxyXdWdJppNmS1wJ/v13tNDNbZOk6tvamHpNMed3yz0t6d7E/it/5Te91X0nvlXRE0nckfbgoekbazq7Ip5NmQr6PVIGVj+cU+3PgnsDbgS8ArwM+D/x4RHxnG9tpZmabcxqTr1u+BfBh4LdHvYmkHwPeA1wM3A+4N/BSYH3cybfzOjZN2H8UePi4Y8zMbPeb5rrliHhDsW9c9fUy4JUR8aLKti9MOr/XijQzs6026rrlkYprmH8cuErSh4rrm/9J0k9Oeq3vx2Zmto90I+Pm/MBm3+YkSZdWnh8qLr8apem65Um+v/h6NvA7wCeAXwTeXVzb/G+jXuhgMzOzWR2OiFOmOXDUdctTKHsUXx0R/6f4/hOSHgz8GvC0US90sJmZ2ZYorlt+POkSr7HXLTcoL/P6bG37Z4E7jXuhx9jMzGzuiuuWy8U4Jl633OA/gCuBH6xtvztpRaqRXLGZme0jaa3Irf3VX1y3fDrwaIrrlotdRyLiSHHMiaTK61bFvrtK+hZwdURcHREh6Y+BsyV9ijTG9ljStP9njDu/g83MzObt6cXX99W2nw2cVXz/c8BrK/teUz8mIl4u6QDwJ8B3AZ8BfnrcxBFwsJmZ2ZxNum65OOY84Lwpjvsj4I9mOb/H2MzMbKG4YjMz20ci2PIxtp3mis3MzBaKg83MzBaKg83MzBbKYne0mpnZgO24jm2nuWIzM7OF4mAzM7OF4mAzM7OF4mAzM7OFstgjiGZmNiBCHM2Xd7oZW8oVm5mZLRQHm5mZLRQHm5mZLRSPsZmZ7SNdxFq3tdPN2FKu2MzMbKE42MzMbKE42MzMbKF4jM3MbB+JEGteBNnMzGzvcLCZmdlCcbCZmdlCcbCZme0jAazlrU09JpH0QEnvkPQNSSHpjNr+20o6T9KVkm6WdKGku1X2nyjpTyVdLumopK9J+jNJ3zXNZ3SwmZnZvB0HXAY8Czha3SFJwNuAuwGPBv4TcAVwkaRbFId9D3AH4LnAPYEnAg8E3jTNyRd7aoyZmW27iLgAuABA0nm13XcD7gecHBH/VhzzNOBq4JeA/x0RlwE/X3nNlyT9DvBOSSdExLfHnd8Vm5mZbacDxdeVckNEdIFV4NQxrzuhOObmSSdwxWZmto90Q6zkm/7Vf5KkSyvPD0XEoSlfeznwVeDFkn4VOAL8FnBH4PZNL5B0K+B/AK+JiM6kEzjYzMxsVocj4pSNvDAi1iX9PPAXwDeBHLgIeBeg+vGSjgPOB75BGnObyMFmZmbbKiI+Dpws6ZbAckRcJ+kSoFoFlqF2QfH0ZyNihSl4jM3MzHZERNxYhNrdgFOAt5f7JB0PXAi0gJ+JiCPTvq8rNjOzfSQQ61Nci7YZRaV11+JpBtxJ0snA9RHxVUm/CBwmTfO/J3Au8LaIeE/x+uOB95AmjDwauEXlUoDrI2Jt3PldsZmZ2bydAnyieBwDnF18/4Ji/+2B15Mmkvwv4A2kqf6le5MuCfhh4AvAVZXH/Sed3BWbmZnNVURcTMNEkMr+/0UKtA29fhJXbGZmtlBcsZmZ7SMRsN7d2jG2nbZtFZuk50n6mKRvS7pO0vmSfrR2jCSdVSyMeVTSxZJ+ZLvaaGZme992dkWeBryKNPD3EKBDWvTyxMoxzwWeDTwTuA9wLfDeYoaMmZnZRNvWFRkRD68+l3Q6cCPwE8D5xYrPvwn8YUT8XXHMk0nh9gTg1dvVVjMz27t2cozteFLFeEPx/C7A7UjXLgAQEUclfZBU5TnYzMw2aTuuY9tpOzkr8lzgk8BHiue3K75eUzvumsq+AZLOlHSppEvXWd2aVpqZ2Z6yIxWbpHNItyc4NSLyjb5PsZr0IYATdGLMqXlmZraHbXvFJullpCvMHxIRX67surr4etvaS25b2WdmZjbWtgabpHPph9rltd1fIQXYwyrHHwQeAHx42xppZmZ72rZ1RUp6JXA6aUHLGySV42ZHIuJIRISklwPPl3Q5aX2w3yPdhO6N29VOM7NFFiHWFnzyyHaOsT29+Pq+2vazgbOK719CWjDzlcCtgUuAn4qI72xHA83MbO/bzuvYJi5oGRFBCrmztro9Zma2mLwIspmZLRQvgmxmto9EQJ4vdk2z2J/OzMz2HQebmZktFAebmZktFI+xmZntMx2PsZmZme0dDjYzM1soDjYzM1soDjYzs30kQuTdbFOPSSQ9UNI7JH1DUkg6o+GYu0t6q6RvSbpZ0r9Kukdl/+0kvUHS1cX+f5P0y9N8Rk8eMTOzeTsOuAx4ffEYIOkuwD8X+x4CfAv4IdKi96XXAycCjwKuA/4r8AZJX4uID447uYPNzMzmKiIuAC4AkHRewyEvAt4TEc+ubPty7Zj7A8+MiEuK538i6TeA+wJjg81dkWZmtm0kZcAjgc9KulDSdZI+JulxtUM/BDxW0ndJyiQ9Cvhu4KJJ53DFZma2jwSQdzZd05wk6dLK80MRcWjK196G1FX5fOD3gf+P1B35V5KORMQ/FMc9FngzcBjoAKvAL0XEJyedwMFmZmazOhwRp2zwtWWqvj0izim+/6SkU4BnAGWwvRA4CXgoKdweDbxe0gMj4t/GncDBZmZm26mswD5b2/454PEAkn4AeCZwciXE/k3SA4rtvzLuBB5jMzOzbRMRa8DHgB+s7bo7cEXx/bHF17x2TM4UueWKzcxsPwmIXFt6CknHAXctnmbAnSSdDFwfEV8FXgK8RdI/Ae8HHkyq1h5dvOZy4EvAqyQ9B/hmse9hpOn/Y7liMzOzeTsF+ETxOAY4u/j+BQAR8TbgTOA5wKdJ3YtPKieORMQ68DOk69fOBz4FPAl4SkScP+nkrtjMzGyuIuJiYGxZGBHnAeeN2f9F4Bc2cn5XbGZmtlBcsZmZ7Schur4fm5mZ2d7hYDMzs4XiYDMzs4XiYDMzs4XiySNmZvtNd2sv0N5prtjMzGyhONjMzGyhONjMzGyheIzNzGw/CYjOYo+xzRRskg4A30Na1PK6iLhuS1plZma2QRO7IiUdL+lpkj4I3Ei6lcBlwNWSvirpNZLus9UNNTMzm8bYYJP028B/AP8NeC/pPjgnk24I9+PAWaSq772SLpR0t61srJmZ2SSTuiLvBzwoIi4bsf9fgP8j6deApwIPAr44x/aZmdk8BWiLbzS608YGW0Q8dpo3iYhV4FVzaZGZmdkmeLq/mZktlKlnRRYzIp8OPBi4DbVQjIj7zrdpZmZms5tluv9rgJ8F3g58FogtaZGZmW0hwX4eY6v5OeBREfGPW9UYMzOzzZpljO1a4PBWNcTMzGweZgm25wMvlnTrrWqMmZnZZs0SbO8BjgWulfQ1SV+uPraofWZmNk/FdWybeUwi6YGS3iHpG5JC0hljjn11ccxzRuyXpHcVxzxmmo84yxjb64EfBl4OXIMnj5iZWbPjSEsvvr54NCqC6r7AlWPe69lAd5aTzxJsDwMeEhGXzHKCKkkPBJ4D3Ju0mPJTIuK8yv7zgCfXXnZJRNxvo+c0M7PtFREXABdA7/f6EEnfB5wLPBR414hj7gM8i5QZ10x7/lm6Ir8KrM5wfJMyxZ8FHB1xzEXA7SuPn9nkOc3MbBeR1AbeBLwwIj434pjjgTcCZ0bEtbO8/ywV228BL5H09Ij40iwnKU2T4sBqRFy9kfc3M7PJNFPHXqOTJF1aeX4oIg7N8PqzgcMR8Wdjjvlz4MKIaKzmxpkl2P4GOAB8XtIq0KnujIgTZj35CKdKuhb4FvCPwO/OmtZmZralDkfEKRt5oaTTgDNId4oZdczpwP8DbOgcswTbMzZyghldCLwV+ApwZ+CFwPsl3btYaHmApDOBMwEOcuw2NM/MzDbpNNIw01VSb4ZlC/gjSb8ZEXcEfpI0WfFI5RiAv5b0kYg4ddwJpg62iHjdDA3fkIh4c+XppyV9HLgCeAQp8OrHHwIOAZygEz1L08xs93sV8Le1be8mjbm9pnj+u8BLa8d8mjT58O2TTjA22CQdHxHfmaqpGzh+koi4UtLXAd/A1MxsHgLU2dq1IiUdB9y1eJoBd5J0MnB9RHyVtJJV9fh14OqI+DxARHwD+EbtGICvRcTE66YnzYr8oqTfk3THMR8gk/TTkt4L/PqkE85C0knAHYCr5vm+Zma2pU4BPlE8jiFNFvkE8ILtOPmkrsgHAC8Cvizp08ClpAvpVoBbk/pA70eauv9i+mVko3EpXjzOAv6OFGR3Bv4nKdn/fraPZWZmOyUiLgamLgsj4s5THDP1+026g/YXgcdK+l7gsaSguy8pgQ+TEvgQcEFETDOB9BTgA5XnZxeP1wFPA+4JPAm4FSncPgA8dp7dm2ZmttimmjwSEV8D/qR4bNgUKf7wzby/mZmNpwDlO92KrTXLyiNmZma7noPNzMwWioPNzMwWioPNzMwWyixLapmZ2QKY5mahe9nUFZuk90v6/Ybtt5b0/vk2y8zMbGNmqdhOA06WdC/gSRFR3k9tGXjQvBtmZma2EbOOsT0UuAvwIUm334L2mJmZbcqswfZ14FTgS8Clku49/yaZmdmWKS7Q3sxjt5sl2AIgIlYi4nGkpbQuJi21ZWZmtivMMsY2MI0mIs6W9BngvLm2yMzMbBNmCbYHk1bg74mIv5X0OTZ4+24zM7N5m+UO2v84YvtngM/MrUVmZral9sI42WZ45REzM1soDjYzM1soDjYzM1soXivSzGw/8Y1GzczM9hYHm5mZLRQHm5mZLRQHm5nZPqKALN/cY+I5pAdKeoekb0gKSWdU9i1J+iNJn5J0k6SrJL1R0p1q73FA0p9KOlwc9w5Jd5zmMzrYzMxs3o4DLgOeBRyt7TsW+M/Ai4qvjwK+F7hQUnVC48uBXwB+CXgAcALwTkmtSSf3rEgzM5uriLgAuABA0nm1fTcCD6tuk/TfSStY3QP4tKRbAk8FnhIR7y2OOR24gnT7tHePO78rNjMz22knFF9vKL7eG1gC3lMeEBFfAz4H3H/Sm7liMzPbZ+ZwHdtJki6tPD8UEYc21BZpGfgT4PyI+Hqx+XZADhyuHX5NsW8sB5uZmc3qcERs+q4uxZjaXwK3An5u060quCvSzMy2XRFqbwLuBfxkRHyzsvtqoAWcVHvZbYt9YznYzMxsW0laAv6aFGoPjoh6WH0cWKcyyaSY6n8P4MOT3t9dkWZm+0mAOlt7CknHAXctnmbAnSSdTLpZ9ZXA3wD3AR4JhKRy3OzGiDgaETdK+gvgJZKuBb4JnAN8Crho0vldsZmZ2bydAnyieBwDnF18/wLgjqRr176HVJldVXk8rvIevwn8Pamy+2fgCPDIiJg49cUVm5mZzVVEXAxozCHj9pXvsQo8s3jMxBWbmZktFFdsZmb7ie/HZmZmtrc42MzMbKE42MzMbKE42MzMbKF48oiZ2T4iIOvudCu2lis2MzNbKA42MzNbKA42MzNbKB5jMzPbT/bBBdoONtsR2dJsf/W661u8HLmZLYxt7YqU9EBJ75D0DUkh6Yzafkk6S9KVko5KuljSj2xnG21rZUvtmUNtM6+b9H6zPsxs99vuMbbjgMuAZwFHG/Y/F3g2aTXn+wDXAu+VdPy2tdC2zDyCYbMBs9OvN7Ott63/QiPiAuACAEnnVfdJEun+O38YEX9XbHsyKdyeALx6O9tq87MVQZAttWfqnpx3G2Y9v9luki34X93dNCvyLsDtgPeUGyLiKPBB4P471SjbuK2ubqZ5/61sg6s3s91pNwVbeWvwa2rbr6nsGyDpTEmXSrp0ndUtbZzNZjt/4Y8613a1weFmtrvs6X+REXEIOARwgk6MHW6OFSb9otfy8lTvE2tr82jOppXtHdced02a7R67KdiuLr7eFvhqZfttK/tsl5tXqDUdO0uwzLMd1dc43Gyv0z64jm03dUV+hRRgDys3SDoIPAD48E41yqYzabxJy8sbCpP6e0xqQ/XrRt9n0mvHvd7jbmY7b1v/BUo6Drhr8TQD7iTpZOD6iPiqpJcDz5d0OfAF4PeAI8Abt7OdNpsNBUm7eE1ntgpnmm7Bmduywfdx9Wa2O233fy1PAT5QeX528XgdcAbwEuAY4JXArYFLgJ+KiO9sbzNtGhuujNrt5u/rxoTeqGCZVDXOpN62WnsmhazDzWxnbGtXZERcHBFqeJxR7I+IOCsibh8RByPiQRFx2Xa20eZjZJU2LshGHT/idZsZr2t8/wnn20g73DVpu6B0xwUAACAASURBVE6A8tjUYxJJx0t6uaQrilWkPizpPrVj7i7prZK+JelmSf8q6R7z+Ij+F2cbMnNlVAkEjXhtTKpu2u3GqmlSt+RQe2YJ1zm1w9Wb7TP/G7gX8GTg68ATgYsk/XBEfEPSXYB/Bl4PPAT4FvBDpKGnTXOw2cxmCrVaiIwKtVH7hsJuxlDZSKiNa2OvPQ1jhA43M5B0DPALwC9ExMXF5rMkPRJ4GmnuxIuA90TEsysv/fK82rCbZkXaHrDRUNNSezAwxnX5Vd+zeF3ja8ede0J7mt5/6Dxj2jPqPaeZNWm24NpAC1ipbT8KnCopAx4JfFbShZKuk/QxSY+bVwMcbDa1Ub+UG3+Zj+p6bAqzaca36u9TO0fZjmnaU77XpBAb166NBG3J4WY7Lcs39wBOKld9Kh5nlu9dTPb7CPB7ku4gqSXpicCPA7cHbkNaEP/5pCUUHwa8CfgrSY+Yx+fzvzCbyrhQGzCu67G6r7q9qXuuKdw6nd77jeoSHGjPiAptqnONUukK1VJ7sKu01k06zfif2R51OCJOGbP/dOD/kMbXcuBfSeF1b/oF1dsj4pzi+09KOgV4BvAPm22cKzbbsA2F2lJ7MNSq2+qPukplNLZ7sqG6mrkrdFybxk2EmbJyc9Vmiywi/j0iHkSqzL43Iu4LLJHG0Q4DHeCztZd9DrjTPM7vf122IVOHWlOVNqIaG1L/5V+t0ipVEzRUcEwOnZHnGWepndpRqRQb2+DKzYyIuAm4SdKtgYcDz42INUkfA36wdvjdgSvmcV4Hm01Ury7mHmrjtlcDr3yPWrBUzxnrncmBNs3Y2iidTj/cymObuianCDfPkrSdoNj6+7FJejipR/By0mpTf1x8/9rikJcAb5H0T8D7gQcDjwcePY/zO9hscyaFWlN3ZKs1+v3y2uqs1fcvg2KKgGtq28juzWm0WqltZWDV27DBcDNbULcE/idwR+B64O+A342IdYCIeFsx4eT5wLnAF4EnRcSmx9fAwWYzGjc5Y5ZQi3b6qk4tyEaFXhkqpaZwaTJLmI0L3HJ/NdzK93e4mQ2IiLcAb5lwzHnAeVtxfgebjTXtJIfG7sfSiFCrf9/4vmXwVUOnGnLVgIMUMtOE2YQQa2qXOvlcwq3K3ZFm8+dgs5HGjq01TdIYN6ZWD7WlMcGy3q/i6gHTCxcYX8U1BWxr+kAdULS1XCFP9XNXx93GXQ5Qvr5WtTncbFsVa0UuMgebNRpbqU1a93FcqJUh0U5XmqjTbXj9iNBZzwcCSdV99ZCrnLfUGGbjArbp2KINvXPXx91q4dYzpmoDh5vZPDnYbEhTqDUuWTXqOrXq81ZrqEorQ63+/TjqdAdDqBJyA1VczVCYjQmySW3ptaEMt6auyVq4TdslCQ43s3lxsNmAiWNq48bQZgi1aQKtWs1Vjx8IuXrA9Y6vBNhSvXKbLkxHHafKeUeOu0FzuJXvMeZ+cg43s81xsFnPLKthjJws0jSlvyHUolUJqryhO5LhYCmDLtpZP/QaAm7AUnVcbfD9pg246vHqdNPXyjkHwg0Gr3Mr216G24Sqzcw2z8Fm81ENxcbZj82hBgw9b6K8OzAuNzRGV3QR9tszukprCrRp2tBrS+2cvXCD6cbbJnDVZlttqy/Q3mkONgOmvB1NZZ3G6vOmmYjVLsimUJtqPKuifN3YgGsYP2sKtEkhNrILstMdrDQrY27A6PG2UZ/R17SZbQkHm83HiC7IUaHWrYRH1jAzsilcqsFSBlxj9ySjK7Rpg7U00E5Gh1vAcJdkqT7WNkUV56rNbOMcbLa5sbVRXZC1sa1RoQYMPW9sY606i1Y2snqrnrf3/ZhAm+b85XFThRu4S9JsBznYbDpNsyHHdUFSBFqtcur2KriBq9B61Gm+cLTbznqVXW8SRxFu1W3Vc/W+bwi1cWE2rm3VcKu+Z/95a/AC7o4nkdgusw8u0Pb92GysxvuJTVjVo9oFCSlY6qHWbfUfVdFW4wMGw6hpvKzpUoJJoTbqXL1jK+2s7qtPgol2NjjGN/DzGL0ayrg7bZvZxjjY9rmpJo2Uz6vHjqrWGrogoxYk9TCrhkd9X/V1MH24Vbc3hVpTiE1qSxlu3UqA1j9fCvVKqE1xF2+Hm9l8uSvSJqv+cm76vjX4y7zsgqwGTred9UJtUjdkU7hleRBt9bsDG2ZNVq+Ha5r52NQNOilI6+3rtkQGdKl0jdbH26h1ScLEGZKjeJzNbHYONms0topomDCStrdGjqvVK7VuJTiyItCqYVIfa+u21Bhu9bE1aK7g6qHWFGiNFVy1nZV2RVsp3BheIaW8eBvo/3w6DctqVRdM9tR/2yaK6P2bW1QOtn1sVDfkxHUhYWwXJAyOqw2OTYluW0QRLMpjIDwgBV09ZHrV0phwG6raJoTapCDrf5bidVBUa6mCLN+7WhtWL94emiUJvaptlmW2zGw2DjYbMLZSa1g9f1QXZH1afX98KoVat51WP4hK5aReWDRXc9OGW9mW9F6TQ21cmKX9RVvQYLhRreCal/8qf04D17bNOCPS3ZFms3Gw2Xj1MFtqN2wb7oKE4XG1aqhFS3Tpd4dMCrmyimsKt0maQq0eZlGfKNIe3tclGsOtOt6W3r+4BAAGqzYYHmvzHbbN5s7BZj1D1Vp9Ca2qslqb0AU5MtTaqQIqjQu5tL853IDhqq1WrUFzqI0Ls/r+sr1jw638OdYWdR64f9sUlZrDzbZUjL5edFE42Gw67XZztcZwF+Q0oQaDQTIu5KCs8AbDLZ27eaZk/Zq53vcNodZUndW3d1v9dk4TbqWBiSSl4j8KAl+sbbYFHGz7VH3iyKhqrVGlWttIqHWH1yoebNuokGMw3MpJHINB13zNXD3UZgmz+vN6uLVG/O93qEuy3NEUYu6SNJsbB5sNGxdqG1QNtXoXIAwGWFU95MowSfuofT8+1HpjZe3pwqx6TLRAOdCwFFEvZBnfJWlm28PBZmOp7H4sv26wWitFS0SL3qzIypl6340LuS4BbZF1+l2A6TWV7ydUafU2lWFWD7Le/t6syM11SQ5UbdXt6w3XuNU/u2dG2hxlW7xWpKT/AL6vYdcFEfEISc8Dfh74QWAV+CjwvIi4bB7nd7DZYDdkpVprnDQywrRdkGWoRatfbaXzDo6nqRiWyluDC7aWFdy46q2//FZzlVbtDp0UZtXtXUB5fzbnuHBjhgwauKbNF23bYrgPUO3Mvz3wceAtxfPTgFcBHyP9/+4FwEWSfjgirt/syR1sNpNR1VrV6FBTP9QyIIOskm7VoKtdFVbsr4+3jajeJlRp1e7QSWEGRVuL90/bm8ONylhb/cLt3qfo1CaSmC2giLiu+lzSU4FvUwRbRDy8tv904EbgJ4DzN3t+B5v1NVVrtW7IuqYuSKAx1PID/VCrVkGlahA0b9dQuKVjU7DkLaE8xlZp9UCbFGYDxwIcgNbqYFvKUOu2NdQl2aov91Wd+l/pWhxVtZntdZIEPBX4y4g4OuKw40n/bG6YxzkdbPvcVMtnVdWqtVK9CxIGw6QXDkWoRcMp8nKCBgxUc73tbcg6qZsyb0G3FWQ5RailLsuMhssKaoHWr+L6524Kst7zalu75TVtUIZb1qle51bsy0esSLJeu1gbxoaYuyNt3jSf69hOknRp5fmhiDg04tiHAXcBXjPm/c4FPgl8ZLMNAwebjVObHTlwOxYGq7XetkoXZPm8Oq6WXjccHv33BBW/5/PiGOX0gq4afllHRJ4CDiDL++eaFGZlkI0NMQaPUSe9bmBJrXzw4u1W0bbqXQAGqrZyDcl8uEuyaf1Is13qcEScMuWxvwp8LCL+rWmnpHOAU4FTI2IuffUONkvGTfGvdUOOq9Zg+Hq16rhalBNIerMLh0U10OiHXS/UaiGXdYrzFiFXDdPq+9XDLEZ0Q9aVbVVZAdKfIFJdGqyVVy4kh6GqbWj9yKYV/0vujrQFIOk2wKOAXx+x/2XA44EHR8SX53VeB9s+NPbmovUbiC7Vq7bhe6z19jV0/zV1QZazI6E2M7KuBVkZbkWINYVcNQizjgaqw+p4XtSCrvc5JlwwXr5m4HW1LkmKsb3qRJKyaitnSPYCbr1hxf8iyMZVbZ7yb3vQGaTp/G+q75B0LvA4UqhdPs+TOtj2mXGhNk7T3bH734+6XqxfrZWh1vu+DKoJoZJX9ve6IOsVXfG7Xtnw5JR6mPWn+U/90Xu6rcrMy7zfJdnK++/b64psq7ciSTlDsnfB9lLLsyNt50QM3ah3KxSTRn4FeHNEHKnteyVwOvBo4AZJtyt2HakfuxEOtn2sN3GkvthxvVob0Q05sLrHuMkaWT/UqtXarOFShmGvWisqumo1F53JYbaRUCuV3YzldXjVSwCqVVsrp7ciSRlw0coG77Bdv7t2vWpr6I501WZ7yGnA3YAnNux7evH1fbXtZwNnbfbEDrZ9ZEPV2sANRQe7Ifvfj67Wyu9hONSiNfvMLOXlDUqL98kZeJ61RgdZ/+v05y3PNyCvjM8VXZL1qq0ca2uVN02tdEf2xtqq42zj2lCbGelws70gIj7A8GI75b7G7fOyq4JN0lnA/1/bfE1E3K7hcNsKI8bWmroh69et9faPqNbKbdGC7vI091EbPEYd9UJJvUAbDrpRQVat1Krvrc7of2PRiv4lCOV1c63yJqPp87V65+xXbVTuPkCe7j7Q6lVr2eA4W6czfJ82M9uwXRVshc+TStiSByPmYFK1Nu7atYG7ZFe6IfvbhpeqGlWtlZVaPbQGzlebKtkLs2K7uv2Kqbw0oBp09SCrnit6MyOjf3wlZNUw9JAmrKTVHlVcUlCt2qI7ONZWXtdWnfqfKrds8jjblLMhXbXZhkVthu4C2o3B1omIq3e6Eftafayt0HgDz9r9zcqZkPVqbSDUWsPhlV5bC7va815g9aq0fiiVgaROP8jqITbuPcn7VWZj92M7LWFcDeqmqq28rm1Ud2TjOFtDmPmaNrON243B9v2SriRNEb0EeP48r2/Yj8ZWa5NuUVNMHClDrT4bEhiYNFJ9Xr1mbSDc2jEyxKrb1VDVRUf9cKoEUC/sssp7VN5rqvdlsOJjIOAE7SjG9NSbPKJK1VZe1zapO3JonK0+iWRKrtrMmu22YLuEdN3D5cBtgN8DPizpRyLim/WDJZ0JnAlwkGO3sZkLqhpyEwKvOhsSmi+GjlrAlV2QkdEYOmpHf6Q5K1YTaQ12mXSzSqnXrQRbeZF25dhqgCmrrDFZf888Q8uR3q8dvfeiFb3qLVoBa1nx+frrUUarX7WpPbo7stfOojtyYJytrlrB+UJts5ntqmCLiHdVn0v6KPBl4MnAOQ3HHwIOAZygE7f2BkN71EavWxtQTBypj681dUOmr2NmQmZld+TkMFMraLUHQyivrHEXlYpqIPDKt6sEmCpBWn3PvJPRauVErt6NSnshx2BQ0ooirNXvdmylQCuvyWvqjiwv1h51p23AAWY2R7sq2Ooi4oikz5CuhbAtMjBxpPp9w2r+pXHdkNWlq6ozE8suyFju9quphjCDFD6tYhXkdrGvk2e0liuhVFm5OG8IjWqAtbLBgOy9Z5al92mn9yjDsinkyIv9RZdked+43ufslp91sDuyvFh7aJytnEBS7Y6ckbsjbSMW/e7uuzrYJB0Efgj4wE63ZdE0rerfpJwR2TRxpNS0Wn516azeo+iCVDvGVmetrEu71aVV7Fsu+vPW8sGgzfN+Wzq1iq1d625s1Z5X37N8nzLkqlUc9EMuWpHa39XA5ypnSHbz0d2R5TjbRJVxNk8gMduYXRVskl5KusncV0ljbL8P3AJ43U62a98ZMb5WvX4N+l2P1e+jFmZpW6VaawXKgqzVHVmdtVrdXvAstXKWspz1boul2oDVei3oqsG3XDu2/lqApSxnqdViPW+xlrdotbrkeVar4jKyohTrlpNAKFYcKWc+tvpjaeO6I9N7DE4giXYLeXkts7naVcEG3JG0WOZJwHXAR4H7RcQVO9qqPWqm8bV6mNWW0aqrj6+lbYOTRtL+4WqtDLUy0Mrqqgy0MoRS8OQcbHVYyYc/y3o2GGxLlTYvZZODrdy+nrV6ry1DLs+zXtdnWcWRBxzowmpG5BoK8XHdkTDieraqMeNs4+7L5u5Is0G7Ktgi4vE73QYbr3phdn18rarfDTlcrakVLC+nX8SjKrQy0JZbOctZh+OWVnvvvdYd/mtb76aE4aqttz3rh8Bat81a3mIlbw9UgGvQq+BSQ9NklW631euSjPZg1aZad+TAz6M9YfJInSeT2BZRhC/QNoP+jEhoDrem6mVUtdZqD4+hjQq0Wy6NupN8cjQvF3Ie3/5jWs3VztF8mbWsxXIr59trB4BaF2VxXG+yCZB3Jv+zKa/fq0/7H1o3spzyv4kJJGY2yMFmfUP3XhszK7JhfK3xuIZqrZV1R3Y7nrC8ynLWYTnLudXSzRzIOhwoKqzVhkrthPbKTB+x7kDWYbXb5mi+zHLW4cj6gYFuz6VW3h+Dy7vkrQwtdyHXyO7I8mJtaAi1aSaQmNmmONj2s0mrjsDA4sf1GZGjxtcGuiEbqrUDy53GKu24pVWWs5xjWmuc0F7hQNbh2NYqBzV9l9xKZYBvmtetRJsb1m/RC7jlLOdI50Bj92TeysjbXSIXeSsb6o6sK8fZet9Xr8GrL63lCSRmc+NgW1BzuTC7Qbf2C7xpfK1pVmRZrR1YSqF27NLaQLfjce1VjmmtcSDrcEJ7pRdot2zdPFP7bjnj5zkQSxxUh+s7t2C12+ZA1uGY1ho3rh/DcivvB1w2vmorfxZN42z1UKuu9N/IK/3bVvIiyLbvNFRx9VmR3fZguA2GWPV1/RVGymoNigkiU4bawWydW2SrbJWbugf4FsdyYvumXvUGcMulo9y4fszAscutvFe15Wv9SSTlSiRNNzAtx9n6z2ecRGJmM3Ow2YY0DHcNhlplIeJylZFybG0pG5wgUoZateuxGmoHlaZwzNIl2aTaTdmTwUp3KV2Y1qU3ngewXLlkYCVvs9TKWWu1YD0FdbctWEsBzlr/Or7qz0OTehiXWjBqzcgZeMq/WZ+DzTasKdygUr2VhV7RDQmDU/DLSSJlqB3IOhxUhwPZei/UbpXdzEGtj23HsVnz/pu7SwPPm97nW134rvYRvtk5DjI4ttWvDtNMyjTr8mAxSxJSQOetDNZJFdt6eUeAykXqlQkkpZmn/JvZhjjYbFPqkya6tS7JaKVFjssp/mU3ZKrW8oFq7cT2TRzI1rlV6+aBUKsG18GmO4H29vVDYyXEwdZwF+ZKw43gBsKtUnGV1dta95jeZJeyO3INesuC9e/7NrJpZrtHxFx6CXYzB5tNpWn1kVK9cquGGvTXgyyvWzvY6jR2QR7I1jmowUrtxNbqQGBVHVS9Tars629diW5le6V7MQRFaFbDbVVL0OpfXlB2SR5stbi5lXOUJVqVcTZVPnN5A9L6z0f58EXbZrY1HGw2F/WltErV8TXoL3VVhkV1XO2gUhfkQa31KrWDiqEAO1bNf20PqN/1uBr9Ku/YIuRujsG+wYOCb3W7kK2zEum1B7N1VmJ9YLyt2iW5lOUpoItxtpxsoBuyrjrl38y2x+j/hpttk5vzA0PbyqBpPD46QyEFKczKR9Nrhs/RrT0ffceDtaKPdb3p+obNWPAuIbOd4GCzQQ3rE6rT3fB1L+XdqMt7p5XBUAZF2d23Em1WYomV7hIrscxKtLm5u8RKiJXoDj2gH3BNQde0b/h9xEpkxXmKdhQTTsrnTaudlLrFOpLKR6+8oo1Ua57daFspQJ18U49ZSHqepJD0isq24yT9qaSvSzoq6fOSfmteH9FdkftZpzN59ZH1HMaMr23G0Xy5t+LHsa1VVrtLHGytc1P3AAdbaynoYp3ekvmFg4qhauugsqFwG67IhgMohVoK05u6B1iJJVa7S9ycH+iFWrkeZbnUVvU+cNPIJuSUVx2xRSXpfsCZwKdqu84BHgqcDnwFeCDwGkmHI+INmz2vKzbrq1UKm/2FW1YyZWVTBsJK3u6t0F+Gx835gTFVW9Z7AEWl1X+kbc1V3fBxWe0xvlor21e2t34PuGlNCjcvgGyLRtItgb8C/htwQ233/YE3RMQHIuI/IuL1pNuU/dg8zu1gs2ab+EU77qLkgVvDVLojV3tBl8IlVU/LrMQSN3f7j3owwXDQDT4Gj6++1829IBus1so21ds58OPpDHZDTrwQG8g6QTZt1+SMt6zxxdm2Cx0C/jYiPtCw70PAIyV9L4Ck+wMnAxfO48TuirRkmm7JMbJuvcMwKX/x550Mivkg5Wr5y1mn1x0JqWo72O6kiSNduEW2OrxaSP0kIy7OrqpfqF2dmLIS7RRqlWOqk1l6t8WhH8qdGbsi67JOF+UN45ZjwmzUTUbNZhfz6CE4SdKlleeHIuJQ+UTSrwJ3BZ444vW/Abwa+KrUW1LomRHxzs02DBxsC6u73tn4Qsi1kFOni9pdsk6XaIvopDtDp3uOBVlH5C1QB7JWyh7loG66OSfdtFjw6no7LanVavXGq45bWuVb68f2VvS/vnOLtPpHC77ZOW5grciDWhtavX8lX6o874dcfVZlPSDLGZBlqN2YH5u6PvMDfLtzsHcrm7Vuq3crm/VucXftbtb7XOSivGZcOWQz/r4Y6u6tVF7hKsx2r8MRcUrTDkk/CLwYODWiYYpy8kxSd+TPAVeQxtheKuk/ImLTVZuDbZ+KtTW0nH65x3oHlSG43unfly3P08XHS6PHlbI83bZFOSgrvuapUlMnbet2RDfLyDvB6lqb5VbOetZ/z+VWzlq3xdF8mVst3dzvllSaPFJWUwcr1Vmq5vrVVD30YHj6/k3dfiXWH09L3Y/l6v6r3TbfWj+WtW5r6O7a63mL1bU2eSdL44arWfE5lb7mDD2yPM2MVB6okx5D8tx3y7ZF8uPAScBnpN6ErRbwQEm/BnwX8D+BX4yI84v9n5J0MvAc5tAd6WCzvmql1tA1mab9Z6hY2Fg50B6eaTjwy71bdEcWVVvezTi6lkJlqZtu4pluMJpeW1Zvq912b6X/VZUhVKnIKl2HB7N1bqLhWrgxXZCrlYkiTVVaGWrfWTvAejeF2tG1pcZqrfys6WeUumWzzuip/r3LJ0Zdw+aQs73tbcCltW2vBb5IquQgDUzU/wHkzGneh4PNhtWqNtZz1M6I1uDfuazDQHdktCCK7sioVm1rGd1W9Kq2TpaxVptd+O21AxxstXrVW7kySXmPNBhcoHhV1XGy0eNsq93B8bRSOY42sUorQm0tb9HJs7HVWlM3ZNYpH2niSFYfV6tWa01djw45m7eILf17FRHfAr5V3SbpJuD6iLiseP6PwB9KOkLqinwQ8CTgufNog4Ntvxs3aaRhX3WcLa8sgKzG7shKJVOr2lbX2uUqVSOV1VsZbNWZitXby1QD76A6jQFWf/1q5XKDepXWFGp5nvWqtVgrQq07/FmbuiHryokj4y6nmGV8zTMibQ96PKk78q+AE0nh9vvAK8a9aFoONgMq42xlmFWqNnVy6LRQu0u0szRWVOmOzPI0eaSs4OpVW7YGXTKCNLGkA3CANF7VylgrKrVq12S9eiu7J4HeRd1N399MvxIrNX1fznacJtA6Rah1Vlt011tFhTZcrdW7IcvsrY6vjZwJWYaTKzRbQBFxWu351cBTtup8DrYFNuvMyMZwA2i1et2RanXJoFe1VeXL6Zf88FST4kLtWriVd6Nu5d2pAg7S4sllKJVh11TR1Z9Xp+2v9Zb1mi7Q8k6q1LrrrdQFuZaRral4QGsNslVorUNrdbBaa60GrZVurxuyN81/Pe9Pua6FWjjkzDbFwbaP9WZGNnVHVrflObRbvWn/QG8SSQYDU/9bqyI/kLoks4HboQ2HG3mkrskxAQdppZKDrSLYWuXdATrpPmm1oCuNCjKgN743daDlWZosMiHUyi7I1loKtOrYWrVaG+iGHBVqtT8nM5ueg80GDEz9h/SLt91GnZxYSuEWraw31tYl3ZomDXOJzoEUbpCmOKVf9OWbVcItF7HcpdvWUMAdoEP5kpuLgbiygksh10kXeLdy1rptlis3BK1aq64eUgmz8v2mDbQouxxHhFov0Fb7XZBZB1qr3V4X5FC1Bv2KrKkym7Ja8/iazSxi4RfadrAtuKm7IysV2lCXZG8cKF3XVq3aMiA60bsou72awi3r9CeWdA8Mh1u0U8CRB9GKwYDrZCwvdwaqt/W81avgyhVAUldlpxday5U7eTYFGfTvLlAGGqQ1LFfX26MDLRfZushW1RszbK1VAm29H2rt1cHr1npdkJVqrTjp0B+BL8g2mw8HmzUaCLfKJJKyaoM0Q7JLv0sSID+QJpOwlm6+GS1gdTjcIgd1RHc5oBVDAbeai/xA3uueBFirdEcutfJexQUMVHNNYVZ+Xw0zSMtjra21RwZaCrLy4uvBUCu7IMvJImWopW7IfhckMFCtDXRDNnVBemzNbFMcbPtcdQWS+lhbL9wqv3R1dA2OScdndCj/CkUuymsrsw7kB0Q3D5SnMbesC+XqV8qhu1yEHql6i1yDAbcMndXKBJOsSyfPaLe6rEGvK7LsolzvtljK8sbKDAYDrVzrsex2zNdajYGmLr0qTTn9UOsMThRprUXR9ZgqtXKyiDpBa6XTn95/dD2FWm2lkXGh5vE1s9k52Gw43GCwW7J2fD3c1MnID7Zp0UUdEe0Ucv1r2FKQqfidrgNp1mC+nKqg7nJxmUAl4AKGxt9aRRdoGXBH84xWq/kGqPWuxvJ15Q1PGyeGVAJNnWql1lylZWvDVZoqlVoZatlKp1+plaFWVGubCTWPr9mGBAvfK+Bgs2aTVvtfz9F6jpZadA+2aa10iHYGB1t0IYVcroHqLV8unhfVmyrhFi0GAm7k+FsRcHk3o5WlUCsDrjehZESgWtmONwAADM5JREFUlWEGNE4MKQMtVWcaCLRJVVoZaq3Vbm9MbWyoNf28zWwuHGwG1Kq2UhFujVUbEO10fVsGdA+mv0qtFdLyW8UFzPnBrH9tVwc6B9KrlRfh1klBNlC9FZXeUPdkJeBgsJuydLQSZkDztH2YOI7WC7VOfyp/fYJIdeZjWaWVoZatdobH1Go/23HXq7kL0mzjHGzWMy7ceqrjbQyGW7QzIs+ANl3SjElW0mUB5Vhbe3Wweuu2+wHXPVDttkwTSxrH39qCLCBn8Dq4SsCNm7YPNI6jlYGW5f1p/NVuxyynsUprHE+rh9qoas2hZjZ3DjYb0BhuVMba2u2hcIPUbV9+X467RTsbqN4gI1qD1VsZNFGuuzjF+FsZcGoH3W5ruJty1LR9SCvyN4yjZWvlaiEp1LLucLfj1FUaDHY/wsC4WvnzHPXzN9tSEQt/aYmDzYYMhduYLsmSKMKtGHcrq7dqwGV5kB/IemNqytPlAb0qrfjbWM6WTBMbi7sG5OMCrhyTYyjQoAjPWqBN6nasVmnlKiKjAq03lga9i6+HZj82hdomZ0BmS21PIDFr4GCzRmPDrby+DQYXSy5fy5iA6wTRVi/gALKW6BAD1RuA2v3qrVy2a3zAFRUc/TG0UYE2bZXWNIW/Gmj1e6uN7Hqcc6iZ2WgOtgU3yyLIdVOHW820AZcfyMg6ZcjVqrdickm1eovW5IDrtWHKQJu2SqvPdpwYaDD1BdgONbP5crDZWOPG3ICR1RsMj7+p0x0IOGj3ZjiOq97Ksbe00PL4gItWjAy0+vR95dNVaa2VfPZAA4ea7U5bfKPR3cDBtqA2U6nVTbU6SfUfSlPAdXKi3RoIOKD3tTq5BES33V+SS+3+yiVqjw647oFAXTUGWn22Y9advkorJ4UMTQyB4ckh5feVn8+4n+tmlX/OHmsz63OwLZh5BtpITauTlOet38utopxgUn6f0aFLm9YKcLAFK93edW/l0fmy0nVypPu8dfN0WUC5rfru2Sq9SSn1UCsnh0wbatXFixurtGqgweRQ2+L/ITvgzPocbAtiWwKtrumOAOX2EYa6J0fcKWBUuGUARYBxYDjcqst4jep6rIdaeZuZ6or89e7HxmvSYLoqbRuvVXPAmTnY9rztCrRRY231cANmqt7ULhdOTncKaBXxVd4GJxsRbl2KbsVi1ZK6+gSRaqj17ptWu+C6OvNx7JJYMPiZZgy17eCAs5EiFn5s18G2R+1IhTbKpK7Jcl/tl2wae0vTHss7BTTeBqcIt/K6tmq4QeqirI67leNpoyaJNIVaVr8+rV6pNa0cUp8IsgtvPeOAs/1oF/12tGntqlCrGtU1Wd1X/wVb3LwUQK1usTRXP7yiJbpEEW7AWuVWOAx2RfYquVqo9e+ZNhhq5Q1By6WxylCrThaZJtRmrdJ24n/LDjjbT7LJh9hukS21d2+oNZjmF351gWDl5Xhbqp7Stv5yVlVZZ7B7UcVbN1VqA6+rhFrWKdtQOV+n8oL1yuLF01ZjU1RqjV2622Qv/f0x26hdGWySni7pK5JWJH1c0gN2uk07bad/IU39y3jaLrhqUFQCpBos5d2nS1knBdMomrCvHo7pPaPyffXczfd5q9uLa+7t9N8l21kRQXe9s6nHJJKeJ+ljkr4t6TpJ50v60THHv1pSSHrOPD7jrgs2SY8DzgVeDPwn4MPAuyTdaUcbtoP28i+irZhM0RRQ42QjTjn1+3SmrNZmsJNVG+y96t/2nNOAVwH3Bx4CdICLJJ1YP1DSY4D7AlfO6+S7LtiA3wbOi4jXRMTnIuKZwFXA03a4Xdtut/zy2Y5fwtUKKRv4vjl8yin7Tdt7r+2OruLqr61Wh9NWa4tgN/z9ssUTEQ+PiNdGxGUR8WngdOC7gZ+oHifp+0iFzBOA9Xmdf1cFm6Rl4N7Ae2q73kNK/n1joX/hlFXPlNXPuMqqaRyt3N6kabxurLzyRnOe6bjTVVtpof+u2W5xPClvbig3SGoDbwJeGBGfm+fJdtvf6JNIs7evqW2/Bnho/WBJZwJnFk9XL+q+5bKtbd42Wu19dxJweOcaQrUt87bzn21r+HPtDt+30w3Yjb7DDe++qPuWkzb5NgclXVp5figiDo05/lzgk8BHKtvOBg5HxJ9tsi1DdluwzaT4QR4CkHRpRJyyw02au0X9XLC4n82fy3aziPgv23k+SecApwKnRkRebDsNOAM4eSvOuau6Ikn/G8yB29a23xa4evubY2ZmGyXpZcAvAQ+JiC9Xdp0G3B64SlJHUodUYf+RpK9v9ry7KtgiYg34OPCw2q6HkWZHmpnZHiDpXPqhdnlt96uAe5EqtvJxJfAy4Cc3e+7d2BV5DvAGSf8C/DPwa8D3AH8+4XXj+nf3skX9XLC4n82fy/Y1Sa8kzYR8NHCDpNsVu45ExJGIuBa4tvaadeDqiPj8ps8fMds1QdtB0tOB55JK1cuA34qID+5sq8zMbBqSRgXL2RFx1ojX/Afwioh46abPvxuDzczMbKN21RibmZnZZu35YFvEdSUlnVWsm1Z97LlZoZIeKOkdkr5RfIYzavtVfNYrJR2VdLGkH9mh5s5kis92XsOf4Ud3qLlTm2aNv73852b7w54OtgVfV/LzpDHG8nHPnW3OhhxHGiN9FnC0Yf9zgWcDzwTuQxpMfq+k47ethRs36bMBXMTgn+HPbE/TNuU0Jq/xt5f/3Gwf2NNjbJIuAT4VEb9a2fZF4G8j4nk717LNkXQW8JiIGLka9l4j6QjwjIg4r3gu0vTeV0TEi4ptx5B+ST4nIl69U22dVf2zFdvOA06KiJ/dqXbNg6TjgBuBR0fE+Yv052aLa89WbPtgXcnvL7p6viLpzZK+f6cbNGd3AW5H5c8vIo4CH2Qx/vz4v+3dT4hVZRjH8e+viVJciGWDMP1xURE2UIvaiJO1SIyCqGjjQssySHAVtIiCgbAMJdoPpYKaFBUY6bgoaGGQZBEZ1SxqnMxIKqgMzYKnxXtuXW63O3PnXjznPff3gbM4Z877nvflnXsfzrnveV5glaTTkqYkTUgaLrtB89Ca428Qxs0yl21go3NeyWX/PT0rH5LSzawFNpH684Gky8tsVJ81xqiO4wcwCawnvWz6BGlZjvckXVpqq7rXmuOv7uNmNVDFF7QHXkQcat4vJh18DWwgvcBuFRcR+5t2P5N0DDgB3A28WU6rutMux59ZDnK+YxuYvJIRcQb4HLiu7Lb0UWOMaj9+ABFxCjhJJmPYIcffQI2b5SnbwDZIeSUlLQBuIC24WhffkL4I/xm/op9j1Gz8ACQtBUbIYAxnyfE3UONmecr9UeR880pWmqQdwNvADDAMPAMsAnaX2a5uFTPqri12LwKulnQz8HNEzEh6CXhK0pfAFPA0cAbYV0qDu9Cpb8U2DrxBCmTLgedJMwffutBt7cYccvxFzuNmAyIist6AzcA0aTnMY8BtZbepD33aT5pSfR74jvQFuaLsds2jH7cD0WbbVfxdpADwPXAOeB8YLbvdvfYNWAgcJgWy86Tf1nYBV5Xd7jn0q12fAhhvOifbcfM2GFvW77GZmZm1yvY3NjMzs3Yc2MzMrFYc2MzMrFYc2MzMrFYc2MzMrFYc2MzMrFYc2Cx7klYXGfSHeqznqKQH+tUuMyuHA5vVwXZga/SeqPdZYJskfy7MMuYPsGVN0kpSHs3X+lDdQdL6Y3f1oS4zK4kDm1WCpPWSfmpdr0zSXkkHOhRdB7wbabHLRplxScclbZA0Lel3STslXSJps6Rvi2u92Hx3VtzxHSQlADazTDmwWVW8Tvp/vLdxQNJi4D7g5Q7lxoCP2hxfXtR1D3A/8CBwALgVWAM8Cmwp6m92FFg9nw6YWTXknt3faiIizkraC2zk38eK64BfgXc6FL2GlDC61RDwcET8AhyXNEkKWCORljz6QtIR4A5SkumGU8CIpIsj4q+eOmVmpfAdm1XJBHCnpCuL/Y3A7lkCzEJShvlWM0VQa/gBmCqCWvOx4ZZyZ0nZ6xd01XIzqwwHNquMiPgU+Bh4SNIocAvwyizFfgSWtDn+Z2v1/3Os9TNwGXAu0qrlZpYhP4q0qpkAngSWAkci4qtZzv8EWNHH64+SgquZZcp3bFY1rwLLgMfpPGmk4TCwqo/XHwMm+1ifmV1gDmxWKRHxG2nyyB/M7d20PcD1km7s9dqSRoCVwM5e6zKz8ngFbascSYeAkxGxaY7nbwOuiIhHerzudmBxRDzWSz1mVi7/xmaVIWkJ6VHgGuCmLoo+B2yRNNRjWq3TwI4eyptZBfiOzSpD0jRpVuLWiHih5OaYWaYc2MzMrFY8ecTMzGrFgc3MzGrFgc3MzGrFgc3MzGrFgc3MzGrFgc3MzGrlb1C1RV5nPOOaAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 576x432 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "from get_var_slice import get_var_slice\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "\n",
    "#------------------------------\n",
    "\n",
    "inputFile = \"case.fds\"\n",
    "varName = \"TEMPERATURE\"\n",
    "pnumber = 2\n",
    "gskip = 1\n",
    "tskip = 1\n",
    "\n",
    "V, X, Y, pdir = get_var_slice(inputFile, varName, pnumber, gskip, tskip)\n",
    "    \n",
    "xlabel = 'y (m)' if pdir=='X' else 'x (m)' if pdir=='Y' else 'x (m)'\n",
    "ylabel = 'z (m)' if pdir=='X' else 'z (m)' if pdir=='Y' else 'y (m)'\n",
    "        \n",
    "#---------- plot results\n",
    "    \n",
    "plt.rc('font', size=14)\n",
    "    \n",
    "plt.figure(figsize=(8,6))\n",
    "plt.contourf(X,Y, V[:,:,-1], 80)\n",
    "plt.colorbar()\n",
    "plt.xlabel(xlabel)\n",
    "plt.ylabel(ylabel)\n",
    "plt.axis('scaled')\n",
    "plt.title(varName)\n",
    "    \n",
    "nx,ny,nt = np.shape(V)\n",
    "Vavg = np.zeros((nx,ny))\n",
    "for it in range(nt):\n",
    "    Vavg += V[:,:,it]\n",
    "Vavg /= nt\n",
    "    \n",
    "plt.figure(figsize=(8,6))\n",
    "plt.contourf(X,Y,Vavg, 80)\n",
    "plt.colorbar()\n",
    "plt.xlabel(xlabel)\n",
    "plt.ylabel(ylabel)\n",
    "plt.axis('scaled')    \n",
    "plt.title(f\"Mean {varName}\");\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
